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\0 This thesis uses Path Integrals and Green's Functions to study 

Gravity, Quantum Field Theory and Statistical Mechanics, 
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■ finite temperature interacting quantum systems in perturbative regime; 

self-interacting fermi models in non-trivial space-time of different dimensions; 
non-linear quantum models at finite temperatures in a background curved space-time; 
3-D topological field models in non-trivial space-time and at finite temperatures; 
thermal quantum systems in a background curved space-time. 
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The sign conventions of the metric of the Part I. is sign{+2) = diag{—, +, +, +) 
with Riemann tensor R'^p^s — ^^r"/?^ — ^^F^/j^ + r^^^r^/j^ — F^^^F^/j^ and Ricci tensor 
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In Part IL metric with signature sign{—2) — diag{+, —, — , — ) is used in order to 
preserve standard formulation theory in real time formalism. 
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SUMMARY 

The problems which are studied in this work belong to three different fields of the- 
oretical physics: gravity, quantum field theory and statistical mechanics. Two well 
known methods of modern quantum field theory, path integrals and Green's func- 
tions, allow one to connect these different branches of physics. It is possible to get 
a number of interesting physical predictions by combining these two methods. This 
combination of path integral and Green's function methods is applied to: finite tem- 
perature quantum systems of different spin in gravitational fields; finite temperature 
interacting quantum systems in perturbative regime; self interacting fermi models 
in non-trivial space-time of different dimensions; non-linear quantum models at finite 
temperatures in a background curved space-time; 3-D topological field models in non- 
trivial space-time and at finite temperatures; and construction of the statistics and 
thermodynamics of thermal quantum systems in a background curved space-time. 
The thesis is divided into four parts. 

The core of the thesis is Part I. It is concerned with the development of ideas of 
local quantum statistics and thermodynamics of ideal thermal quantum systems. The 
goal is to apply the methods of quantum field theory to thermal quantum systems, 
and so to extend and to improve them as convenient mathematical tools for describing 
the thermal behavior of quantum systems of different spins in external gravitational 
field. Particular attention is devoted to the development of the connection between 
Green's functions of the models and their distribution functions. For the description 
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of thermal quantum systems in an external curved space-time, Schwinger-Dewitt and 
momentum-space representations of Green's functions are introduced. The problem of 
introducing temperature into field models in curved space-time for these two methods 
is discussed. The equivalence and difference between these two formalisms are stud- 
ied. The aim of introducing these different representations is to show how to introduce 
chemical potential in the momentum-space formulation of partition functions of quan- 
tum systems and to study their thermal behavior in external gravitational fields. The 
concept of local thermodynamics is introduced, and densities of Helmholtz free ener- 
gies and grand thermodynamical potentials for bose and fermi gases are computed. 
It is shown how average occupation numbers of bosonic and fermionic quantum sys- 
tem depend on external gravitational fields. Low temperature behavior of bose and 
fermi gases is explored, and the phenomenon of Bose condensation in thermal bose 
systems in a background curved space-time is studied. Thermodynamical properties 
of a photon gas in a background gravitational field are considered. 

In part II finite temperature interacting quantum fields are explored from the point 
of view of thermofield dynamics. Two loop renormalizations of a self-interacting A(/7^ 
scalar field at zero and finite temperatures are discussed and analyzed. On this basis, 
the two loop renormalized finite temperature Green's function of a boson and non- 
relativistic Hamiltonian of a boson in a heat bath in an external gravitational field are 
computed. The interesting problem of non-equality between inertial and gravitational 
masses of a boson is studied in this part. 

In part III non-linear spinor models in 3-D and 4-D space-time are considered. 
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The Gross- Neveu model and the 3-D Heisenberg model in non-trivial space-times with 
different topologies are studied. The problem of the dynamical mass of composite 
fermions is discussed. The 4-D Heisenberg model in a background curved space-time 
at finite temperature is explored. The influence of curvature on the dynamical mass 
generation of fermions is considered. 

Part IV is devoted to the analysis of topological field models in non-Euchdean 3-D 
space-time and finite temperature effects. The properties of vector and tensor field 
models with Chern-Simons mass term are considered. The effective Chern-Simons 
terms of vector and tensor fields appearing as the result of interaction of 3-D fermions 
with external gauge vector and tensor fields are discussed. The ideas and methods 
developed in part I are applied to the analysis of topological models with non-trivial 
space-time structure. The influences of the topology of 3-D space-time and of thermal 
behavior of interacting fermions on generation of the effective Chern-Simons terms 
are studied and analyzed. 

New results contained in thesis 

Finite temperature field theory in fiat space-time has been studied in a number 
of journal publications (references in the text of part I) and books [ Bonch-Bruevich 
et al 1962, Abrikosov et al. 1963, Popov 1987, Kapusta 1989 et al.]. Quantum field 
theory in curved space-time has been developed by many authors (journal references 
in the text of part I) and is well described in the books [De-Witt 1965, Birrell & De- 
vies 1982, Fulling 1989 et al.] Finite temperature field theory in curved space-time in 
the Schwinger-DeWitt formalism [De-Witt 1965,1975] is considered in many journal 
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publications (c.f. references in the text of part I). Most important are the works of 
Dowker, Kennedy, Denardo et al. [Dowker & Kritchley 1977, Dowker & Kennedy 
1978, Denardo & Spalucci 1983]. The momentum-space method in applications to 
a field theory in curved space-time was discussed in the works of Bunch, Parker, 
Panangaden [Bunch & Parker 1979, Panangaden 1981]. This method was used in the 
study of the A</7^-model and the quantum electrodynamics in perturbative regimes 
in curved space-time. In Part I, on the basis of the above mentioned works, the 
author develops the theory of local thermodynamics and a statistics in curved space- 
time. The author shows that the thermodynamical potentials (free energy, grand 
thermodynamical potential) of thermal quantum systems may be computed directly 
through the thermal Green's functions of corresponding quantum fields in chapter 
III (III.24) and (III.52); chapter VII (VII.22) and chapter IX (IX. 15). In chapters 
VII and IX the author shows that the densities of free energy of thermal bose and 
and fermi systems in curved space-time can be derived with the Green's function 
in the Schwinger-DeWitt formalism. The author shows the equivalence between the 
Schwinger-DeWitt and momentum space methods in applications to thermal quan- 
tum systems in chapter VII (VII. 25) and chapter IX (IX. 18). The author analyzes the 
difficulties of the computation of the grand thermodynamical potentials for bose and 
fermi systems using the Schwinger-DeWitt formalism (chapter IX) and computes the 
grand thermodynamical potentials of ideal bose and fermi systems in gravitational 
fields with momentum-space method in chapter XII (XII. 1 and XII. 2) and chapter 
XIII. The author considers low temperature properties of these systems in chapters 
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XII and XIII, and studies the behavior of the chemical potentials of the bose and 
fermi systems with respect to space-time curvature (XII. 19) and (XIII. 20). The au- 
thor analyzes the phenomenon of Bose condensation (section XII. 3) and finds the 
variation of the critical temperature of condensation with respect to the variation 
of the curvature (XII. 25). The author studies the thermodynamical properties of a 
photon gas in curved space-time with momentum space-methods in chapter X. The 
author computes Bose-Einstein and Fermi-Dirac distributions in curved space-time 
in chapters XII (XII.14) and XIII (XIII.IO) 

The temperature properties of self- interacting and gauge fields are studied in imag- 
inary and in real time formalisms in the works of Kislinger , Kapusta, Donoghue et al., 
from the point of view of their renormalizability [Kislinger et al., 1976] and finite tem- 
perature behavior of the constant of the interaction [Kapusta 1979, Fujimoto et al., 
1987]. In the work of Donoghue the gravitational field was also introduced in the elec- 
trodynamics with the Tolman's shift of temperature, and the non-equality between 
inertial and gravitational masses of fermion was also discussed [Donoghue et al.,1984]. 
In part II of this work the author computes the finite temperature effective mass and 
renormalized Green's function of boson in a heat bath in two loop approximation 
with a real time formalism in chapter XV (XV. 21), finds the non-relativistic finite 
temperature Hamiltonian of bosons in a weak gravitational field (section XVI. 1), and 
gets the ratio between the inertial and the gravitational masses of a boson (section 
XVI.2) (XVI.9). 

Twisted fields were introduced by Isham in 1978 and field models in space-time 
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with a non-Minkowskian topology were studied in the works of Ford, Toms et al., [Ford 
1980, Toms 1980]. Non-hnear 3-D field models with 75 symmetry and a large number 
of flavors and non-linear models in 3-D and 4-D space-time were studied in 

works Gross & Neveu 1974, Bender 1977, Tamvakis, Kawati & Miyata 1981 et al. In 
part III the author applies the idea of twisted and untwisted spinor fields to obtain the 
dependence of the dynamical fermionic mass with respect to a non-trivial topology 
of space-time in chapter XVIII (section XVIII. 2). The author analizes the behavior 
of the dynamic fermionic mass with respect to the curvature of space-time and the 
thermal behavior of interacting fields (section XVIII. 4). 

The idea of computing the induced Chern-Simons action for vector and tensor 
fields as the result of the interaction of 3-D fermi fields with external vector and 
tensor fields was studied in the works of Redlich, Das, Ojima et al. Das proposed 
the method of the computation of finite temperature induced Chern-Simons action 
of vector type with the method of derivative expansion [Babu & Das 1987]. Ojima 
computed the Chern-Simons action of vector and tensor fields using the path integral 
method [Ojima 1989]. In Part IV the author, using the above mentioned works 
computes the induced Chern-Simons mass term of a vector field in space-times with 
topologies X and x mobius strip in section (XX.2). In chapter XXI the 
author computes the finite temperature gravitational induced Chern-Simons mass 
term using the momentum space method developed in part I. 

In the list of references the author includes only the publications most important 
for undestanding the text. 
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PART I 



LOCAL QUANTUM STATISTICS 
AND THERMODYNAMICS 
IN CURVED SPACE-TIME 

Introduction 

The statistics and thermodynamics of quantum systems in gravitational fields have 
the attracted attention of physicists for a long time. The equilibrium distribution of 
a relativistic quantum gas in a gravitational field has been studied using the kinetic 
Boltzmann equation, and it was found that the solution of a functional Boltzmann 
equation is a relativistic Maxwell distribution for a certain type of gravitational field 
Uhernikov 1962| ]. Later, different variants of non-equilibrium statistical mechanics 



in classical general relativity as a generalization of the standard statistical equations 
for classical systems in curved space-time were proposed [[Israel fc Kandrup 1984| , 



Ignat'ev 1985 



At the same time field-theoretical methods were developed and applied to sta- 
tistical mechanics and thermodynamics. Since the middle of the 50 's, significant 
progress was made in quantum theory of many body systems, and was connected 
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with the development of finite temperature quantum field methods in statistical 
physics. Feynman [[I'eynman R.1953 | studied the A-transition in Helium using the 



partition function in the form of a path integral in quantum mechanics. Martin and 
Schwinger [[Martin fc Schwinger 1959(| considered a many-particle system in the con- 



text of quantum field theory in order to treat multiparticle system from the quantum 
field-theoretical point of view. They described the microscopic behavior of a multipar- 
ticle system using Green's functions, and found that the finite temperature Green's 
functions are related to intensive macroscopic variables when the energy and number 
of particles are large. Matsubara [[Matsubara 1955(] proposed finite temperature per- 



turbation techniques similar to perturbative quantum field theory and calculated the 
grand partition function by introducing Green's functions in imaginary time formal- 
ism. The works of Feynman, Matsubara, and Martin & Schwinger created the basis 
for understanding the close connection between Euclidean field theory and statistical 
mechanics. 

In the 60's convenient methods using finite temperature calculations for interacting 
thermal systems with a propagator formalism in the perturbative regime were pro- 
posed in Ponch-Bruevich fc Tyablikov 1962| , [[Abrikosov et al. 1963|| , [[Fradkin 1965 



and [[Symanzik 1966||. These works stressed the analogy between Euclidean Green 



functions and distribution functions in statistical mechanics. Analogies between sta- 
tistical physics and quantum field theory at finite temperatures were further strength- 
ened in the study of infinite equilibrium systems of scalar and spinor fields 
[polan fc Jackiw 1974||. These works provided an opportunity to construct the sta- 



tistical mechanics and thermodynamics of quantum systems in the language of finite 



temperature field theory, and pointed the way towards a construction of a statistical 
mechanics and thermodynamics in curved space-time. 

In the 70's, interest in quantum statistical processes in general relativity was 
aroused by the works of Hawking [[Hawking 1975|| , and Bekenstein [[Bekenstein 1973 



on the thermodynamics of black holes, where an intimate connection between ther- 
modynamics and the structure of space-time was pointed out. 

Parallel to the research on the thermodynamics of black holes, quantum effects 
present in the early Universe were thoroughly discussed by Guth et al. ||Guth 19"66| , 
Starobinsky 198^ , [Linde 1984|| . From these works one may conclude that, for early 
time and high temperatures, the quantum statistical properties of matter and strong 
curvature of space-time become significant. 

With the development of a quantum field theory using functional integral methods 



at finite temperature [[Weinberg 1974| , [ Faddeev fc Slavnov 1991| , [[Kapusta 1989 



Popov 1987[[ and of a formalism of a quantum theory of field in curved space-time 



DeWitt 1965[ , [Fulling 1989|[ , new possibilities for studing the behaviour of thermal 



systems with non-trivial geometry and topology were discussed. This approach was 
developed by Denardo et al. [Penardo fc Spalucci 1983| , [[Nakazava fc I'ukuyama 1985 



Hu [Pu 1981"^ introduced time dependent temperature T(t) = a{t)/T for an 
equilibrium gas of scalar particles which is described by a conformal scalar field in a 
closed Robert son- Walker (RW) Universe in a different way, by using the equilibrium 
temperature of a flat space-time T. 
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The problems of a thermodynamical description of a non-equilibrium thermal 
quantum gas and conditions for thermal equilibrium in an expanding Universe were 
studied by Hu |[Hu 1982[]. However, in general this problem has not been solved 



satisfactorily. Difficulties in the construction of a thermodynamics of thermal sys- 
tems in dynamical space-time are connected with such problems as the definitions of 
temperature, energy spectrum, and the vacuum state of a thermal system. 

In Part I of this work we extend the results of finite temperature field theory 
in order to construct a statistical mechanics and thermodynamics of bosonic and 
fermionic quantum systems in an external curved space-time. 

The methods developed here are based on using the language of a finite temper- 
ature Green's function in an arbitrary curved space time for the definition of the 
thermodynamical potentials of thermal quantum systems. 

Green's functions of matter fields in curved space-time are nonlocal objects that 



are well defined in a small range of the space-time manifold IpeWitt 1965 



Nakazava &: Fukuyama 1985|| . Thermal properties of quantum systems may be con- 



sidered the properties of the set of quasi-equilibrium sub-systems, which are elements 



of the whole quantum system fKulikov fc Pronin 1993|| . This suggestion will lead 



us to the definition of distribution functions of bose and fermi systems, and some 
interesting thermodynamical consequences. 

The division of the entire system into quasi-equilibrium sub-systems may be done 
in a simple way using Riemann normal coordinates [[Fetrov 1969[| . These coordinates 



permit us to rewrite Green's functions in a momentum space representation, and to 
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introduce a local temperature in order to write thermodynamical potentials through 
these Green's functions. As a result, the thermodynamical potentials will be written 
as a series expansion in powers of the curvature tensor at a selected point of the 
curved space-time manifold. Therefore, at any point in space-time the coefficients of 
the series will change because the curvature will change for all other points on the 
manifold. 

Therefore, in accordance with our suggestions, densities of thermodynamic po- 
tentials will be functionals of curvature, and they will be directly connected to the 
temperature Green's functions by a simple relation. The aim of this part will be 
to develop the mathematical formalism of temperature Green's functions in exter- 
nal curved space-time and to use the formalism to construct a quantum statistics 
and thermodynamics of thermal bosonic and fermionic ideal quantum systems and a 
thermodynamics of a photon gas in an arbitrary curved space-time. 

Part I is organized in the following way: A short review of statistical mechan- 
ics is presented in chapter I. In chapter II functional integral methods are applied 
to non-relativistic and relativistic many-body systems in Euclidean space-time. The 
purpose will be to generalize these methods for a description of statistical systems in 
curved space-time. Finite temperature Green's functions are introduced in chapter 
III for computation of the thermodynamical potentials of quantum gases. Finite tem- 
perature gauge fields are studied in chapter IV. An introduction to quantum fields 
in curved space-time is presented in chapter V. The bosons and fermions in external 
gravitational fields are studied in chapters VI and VIII. The thermodynamics of bose 
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and fermi gases in curved space-time is studied in chapters VII and IX. The thermo- 
dynamics of vector bosons is considered in chapter X. Renormahzation problems are 
considered in chapter XL The phenomenon of Bose- Einstein condensation and the 
low temperature properties of fermi gas are studied in chapters XII and XIII. 
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Chapter 1 



QUANTUM FIELD METHODS 



IN STATISTICAL PHYSICS 



The customary approach to many-body theory is the method of second quantization 
|Abrikosov et al. 1963| , | Fetter fc Walecka 1971|| . This chapter is a short introduction 
to statistical mechanics of simple quantum systems using this method. 



1.1 Equilibrium statistical mechanics 

Statistical mechanics deals with three types of ensembles: 

1) The microcanonical ensemble is used to describe an isolated system which has 
a fixed energy E, particle number A^, and volume V. 

2) The canonical ensemble is used to describe a system in contact with a heat 
reservoir at temperature T. The system can freely exchange energy with the reservoir 
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and T, N, and V are fixed variables 

3) The grand canonical ensemble is used to describe a system which can exchange 
particles as well as energy with reservoir. In this ensemble T, V and the chemical 
potential /i are fixed variables. 

The main aim of statistical mechanics is to derive the thermodynamic properties 
of macroscopic bodies starting from the description of the motion of the microscopic 
components (atoms, electrons, etc.). 

To solve the problem it is necessary to find the probability distribution of the 
microscopic components in thermal equilibrium (after a sufficiently long time), and 
deduce the macroscopic properties of the system from this microscopic probability 
distribution. 

Following this scheme let us consider a classical Hamilton system with 2N degrees 
of freedom in a box of volume V. The classical equations of motion are: 



where H is the Hamiltonian, qi and pi {i — 1, 2, N) are the set of coordinates and 
momenta of the system. 

Let A{q,p) be an arbitrary measurable quantity (an observable). The equilibrium 
average of this quantity A is defined by 



dH 



dH 
dqi 





Pi = 




o 
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= J dqdpA{q,p)g{q,p) (1.2) 

where g{q,p) is the equihbrium probabihty. 

This function is never negative and satisfies the normahzation condition 

J dqdp g{q,p) = 1. (1.3) 

The fundamental hypothesis of equihbrium statistical mechanics is that g follows the 
canonical distribution 



g{q,p) = Z-'exp[-(3H{q,p)], (1.4) 



where 



Z ^ j dqdpexp[-f3H{q,p)] (1.5) 

is the partition function, and = T is the absolute equilibrium temperature. 

It is useful to find a characterization of the canonical distribution that distin- 
guishes it from all other possible probability distributions. It is convenient to intro- 
duce the entropy of a distribution g defined as follows: 

S[g] — — j dqdpg{q,p)hig{q,p) — — < lag > (1.6) 
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Entropy has the following properties: the more ordered is the system the smaller is 
the entropy (i.e., the more concentrated the probability distribution in a restricted 
region of phase space); the more disordered is the system (i.e., the more uniform the 
probability distribution), the larger is the entropy. 

For relativistic quantum statistical systems the equation ( |1.4|) will be turned into 
the equation for statistical operator 

g = Z-^exp[-(3{H - jjN] (1.7) 

where H is the Hamiltonian and N is the number operator. This operator is Hermitian 
and commutes with Hamiltonian H. The parameter fi is the so called chemical 
potential. 

The ensemble average of an operator A will be 

<A>= Z-'^TT[Ag] (1.8) 

The factor Z will turn into a so called grand canonical partition function of the form 

Z = Tt exp[-l3{H - ^N] (1.9) 

This function is the most important in thermodynamics. 

The average value of the energy U may be written from ( |1.8|) as 
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U=<H>= TT[Hg] 



(1.10) 



and the entropy as 



5" = — Tr[^ln^] 



1.11) 



From the equations (|1.9| ) and (|1.11|) we get a useful relation 



S = — Tr[^ In^] 
= -TT[g{-\nZ - (3H + 12N)] 

or 

{l/f3)\nZ = U - S/f3 - fiN. (1.12) 

The quantity 

n = -(l//5)lnZ (1.13) 

is the grand thermodynamical potential of the grand canonical ensemble. 

For canonical ensemble (// = 0) we can write the equivalent equation for thermo- 
dynamic potential F (free energy). 
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It is easy to find that 



F -Vt^ HN. (1.14) 

The grand partition function Z = Z{V,T,iJ,) is the most important function in ther- 
modynamics. All other standard termodynamic properties may be determined from 
this function. For example, the pressure, particle number, entropy, and energy (in 
the infinite volume limit) are 



„ dlnZ rr^dlnZ 
P = T—, N = T— 

S = T^-^^, E=-PV + TS + jiN (1.15) 
oT 



1.2 Statistical mechanics of simple systems 
Formalism of second quantization 

Now we can apply the formalism of statistical mechanics developed above to non- 
interacting many body quantum systems in the frame of the method of second quan- 
tization. 

1. Bosonic quantum system at finite temperature 

Let us study a bosonic quantum system. Each quantum state e of the system may 
be occupied by bosons. Let n (n = 0, 1, 2, 3, ...) be number of bosons in this state, and 
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\n > be the wave function of this state. We will call a state |0 > a vacuum state. One 
may introduce creation a and annihilation a"*" operators with commutation relation 

[a, a~^] = aa^ — a^a = 1 (1-16) 

The action of these operators on eigenstates is 

a\n >= ^/n\n — 1 >, a'^\n >= + l\n + 1 > (1-17) 

and on vacuum state is 

a|0>=0 (1.18) 



From ( |1.16| ), ( |1.17| ) one can get the number operator 



a ' a: 



N = a'^a\n >= n\n > (1-19) 

We may build all states \n > from vacuum |0 > by repeated applications of the 
creation operator 

\n >= (n!)-i/2^a+)"|0 > (1.20) 

The Hamiltonian of the system in the state with energy e may be constructed as a 
product e with number operator (up to an additive constant) in the form 
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H = e{N + 1/2) = e(a+a + 1/2) = ^e(a+a + aa+) (1.21) 



The additive term e/2 in (|1.21| ) is the zero-point energy. Usually this term may be 



ignored and we get H = eN . 

The grand partition function according to ( |1.9| ) will be 



Z = Trexp[-/3(if - iiN)] = Trexp[-/3(e - ^l)N] 



^ < n| exp[— /3(e — iJt)N] \n >= ^ exp[— /5(e — fi)n] 

n=0 n=0 



;i-exp[-/3(e-^)])-i (1.22) 



Inserting ( |1.22| ) into ( |1 . 1 5| ) we get the average number of bosons in the system with 
energy e 

iV=(exp[/3(e-/i)]-l)-i (1.23) 

and the average energy in the form e = eA^. 

2. Fermionic quantum system at finite temperature 

For a fermionic quantum system there are only two states |0 > and |1 > with 
energy e. 

The action of the fermion creation and annihilation operators on these states is 



30 



q;+|0 >= |1 >, a\l >= |0 > 



and 



q;+|1 >= |0 >, a\0 >= |1 > 



Thus, these operators have the property that their square is zero 



a a — aa — 



The number operator is TV = a'^a. One has 



Ar|0 a+a|0 >= 0; Ar|l >= >= |1 > 



Creation and annihilation operators satisfy the anticommutation relation 



{a, a~^} = aa^ + aa^ = 1 



The Hamiltonian for this system may be taken in the form 



H = - aa+) = e{N - 1/2) 



The partition function for the fermionic system will be 
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1 

Z = Tiexp[-P{H - iiN)] = J2<n\ exp[-/3(e - ii)N]\n > 

n=0 

= 1 + exp[-/3(e - //)] (1.29) 

The average number of fermions is 

N = (exp[/5(e - /x)] + 1)"^ (1.30) 

The average energy of the system in the state with energy e is e = eN. 

The next section describes an alternative approach, the method of functional 
integrals for studying the behavior of statistical systems. 
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Chapter 2 



PATH INTEGRALS 



IN STATISTICAL PHYSICS 



Functional integration, introduced several decades ago Ijl'eynnian fc Hibbs 1965||, is 



one of the most powerful methods of modern theoretical physics. The functional 
integration approach to systems with an infinite number of degrees of freedom turns 
out to be very suitable for the introduction and formulation of perturbation theory 
in quantum field theory and statistical physics. This approach is simpler than the 
operator method. 

Using functional integrals in statistical physics allows one to derive numerous in- 
teresting results more quickly than other methods. Theories of phase transitions of 
the second kind, superfluidity, superconductivity, plasma, and the Ising model are ex- 
amples of the problems for which the functional integration method appears to be very 
useful. If an exact solution exists, the functional integration method gives a simple 
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way to obtain it. In other cases, when exact knowledge is unobtainable, the appli- 
cation of functional integrals helps to build a qualitative picture of the phenomenon, 
and to develop an approximate solution scheme. 

The functional integration method is suitable for obtaining the diagram techniques 
of the perturbation theory, and also for modifying the perturbative scheme if such 
a modification is necessary. The extension of functional integral techniques to back- 
ground curved space-time allows one to take into account the gravitational field by 
considering the statistical and thermodynamical properties of the systems. 

The aim of this chapter is to introduce the finite temperature functional integral 
approach to statistical mechanics and local thermodynamics in curved space-time. 

2.1 Partition function in path integral formalism 

This paragraph discusses the scalar field which is described by the Schrodinger field 
operator (fi{x) where x is the spatial coordinate. We denote eigenstates of (fi{x) by 



where ip{x) is a c-number function. Let 7t{x) be its conjugate momentum operator. 
The usual completeness and orthogonality conditions may be written as: 



if > 





j dip\ip) {(p\ = 1 
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{Va I [Pa{x) - </?6(f)] (2.2) 

where 5-is the Dirac delta function. 

Similarly, the eigenstates of the conjugate momentum field operator are labeled 
by Itt > and satisfy the equation 

7r(f)|7r >= 7r(f)|7r > (2.3) 

with eigenvalies 7r(x). The completeness and orthogonality conditions are 

j dn Itt) (7r| — 1 

(7rj7r6) = 5[7r„(f)-7r6(f)] (2.4) 

the overlap between an eigenstate of field operator and an eigenstate of the momentum 
operator is 

{ip \ n) = exp 

For dynamics, one requires a Hamiltonian expressed as a functional of the field and 
its conjugate momentum: 

H = [ £xH{7r, (p) (2.6) 



(2.5) 
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Suppose that a system is in state > at a time t — 0. After time tf it will evolve 
into exp{—iHtf)\ipa >■ The transition amplitude for going from state \ipa > to some 
other state \(fib > after time tf is thus 



< (pb\ e:icp{-iHtf)\(pa > 



(2.7) 



For statistical mechanics purposes, consider the case in which the system returns to 
its original state after time tf. Divide the time interval [0, tf] into N equal steps of 
duration At ^tf/N. Then 



{iPa\ex.p{-iHtf) \ipa) = jlim^y H y^"^'^ 



{fa I ttat) (TTivl eiq){-iHAtf) \(fiN) {(fiN I ttat-i) (ttat-iI exp(-ii7Ai/) \(fiN-i) x ... 



X {(fi2 I TTi) (tTiI exp(-iifAi/) {(pi \ (pa) 



(2.8) 



It is known, that 



(2.9) 



and 



< (pi+i\ni >= exp 



z/d3...(f)^.,,(f) 



(2.10) 
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For At ^ one can expand 



{7ii\exp{-iHAt) \ipi) ^ (TTil (1 - iHAt) \ipi) 



(vTi I Lfi) (1 - iHiAt) = (1 - iHiAt) exp 



i / (i'^x7rj(a;)v9j(a;) 



(2.11) 



where = H{7ii, ipi). Inserting (U), (|2l0| ) and (|TT]) into gives 



((/?a|exp(-zi7t/) \(fa) = lim / JT I d(fi^]5{(pi - if a) 

\ ^ f 1 

exp < —lAt ^ / c/^x [H^-Kj, ipj) — 7ij{ipj^i — ipj)At] > 



X 



(2.12) 



where ipN+i = <fi = ^a- 



Taking the continuum hmit of ( p.l2| ) , gives the result 



ip{x,tf) = Lpa(x) 



exp < 



i J dt J d^x 





vr(f , t)^^^ - H{7r{x, t), if{x, t)) 



(2.13) 



The integration over iT{x,t) is unrestricted, and integration over ip{x,t) is such that 
the field starts at <fa{x) at t = and ends at y9a(x) at t = t/. Note that all operators 
are gone and r one can work only with classical variables [[Kapusta 19891). 
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2.2 Partition function for bosons. 



The partition function in statistical mechanics is expressed as 



= Trexp{-/3 {H - ^N)] 



j dcpa {(pa\ exp {-f3{H - IJ,N)} \(pa) 



(2.14) 



where the sum runs over all states. 

The aim now is to express Zjj in terms of a functional integral. 

First introduce the imaginary time r = it. The limits of integration on r are 
[0, P], then —itf = P and 



Zf3 = J Dn J D(px 



periodic 



exp \ J '^^ J '^^^ 



(2.15) 



Integration over the field is constrained so that 



ip{x, 0) = ip{x, (3) 



(2.16) 



This is a consequence of the trace operation. 
In the equation (|2.15| ) make the replacement 
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i/(7r,^)^/f(7r,(^)-/iiV(7r,¥;) (2.17) 

where A^(7r, (f) is conserved charge density, if the system admits some conserved 
charge. 

Rewrite the equation (|2.15|) in a more convenient form. 



For this purpose, consider a scalar field, which is described by the Lagrangian of 
the form 



L = -i(a^)2-lmV'-f/(^) (2.18) 



where U {ip) is a potential function. 

The momentum conjugate to the field is 



d{do(p) 

Here the Hamiltonian may be written as 

H = 7rdQ(p — L 



dL 

TT = = 9ov? (2.19) 



^n' + ^{V^f + ^mV + U{^) (2.20) 



For evaluating the partition function, return to the discretized version of ( |2.15| ): 
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Using the equation for a Gaussian integral 



1 °° 

/ dx exp < —ax^ + ihx 



exp 



1^ 
2 a 



one can write 



f dwj ( At r 

j exp ynij[Vj+i - Vj) - —tt; 



1 / i^j+i -Vjf 

exp < — - 



and after momentum integrations in ( p.21| ) one gets 



Zp = N' jirn^ j (\{ d^^i exp | Ar | 



1 ( - 

2 I Ar 



Returning to the continuum hmit, one obtains: 



/ /3 



Zp = N J Dip exp 

periodic 



dr J d^xL{ip, dip) 
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Equation ( |2.25D expresses as a functional integral in time interval [0,/5]. The 
normalization constant iV' is irrelevant, since it does not change the thermodynamics. 

Now it is seen how this method works in application to a non-interacting scalar 
field. Let the Lagrangian of the scalar field be 



L{^M) = -\{d^f-\m\'' (2.26) 



Define the finite temperature action of the Bose system at temperatute T in a large 
cubic volume V = as 

5^ = -- 1 rfr 1 d^'x^ (dl + 92 + m^)^ (2.27) 
V 



Due to the constraint of periodicity of the field (p{x,T) (x G V, r G [0,/3]) (|2.16| ), it 
can be expanded in a Fourier series as 

00 

ifix, t) = (/3\/)-i/2 J2 E ^nik) exp likx + cu^t) (2.28) 

Jl= — CO ^ 

where (pn{k) is the Fourier coefficient, Un = 27Tn/(3, and k = 27Tm/L n, m are integer 
numbers. 

Using the identity 

J dTexpi{uJn - UJm)r = f35n,m (2.29) 
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one obtains the equation for the action in terms of the Fourier coefficients 



-| oo 

^/^ = -2 E j:v:ik)U + k' + m')Mk) (2.30) 



where ip_n{—k) = '^*n{k) goes from the reahty of the field (p. 
Then the partition function (p.25| ) may be written as 



Zp = N' j WWd^nik) exp {-iy.:(A;)(u;2 ^ p ^ m')d^n{k)^ (2.31) 



where the exphcit form of the measure of the functional integration ( |2.31| ) is 



D^ = \{\{d^n{k) (2.32) 



So, integration over the coefficients ^n{k) in ( p.31| ) gives the following equation for 



the logarithm of the partition function [Bernard 1974] 



In^, = -\\riDet{u^l + ef) = ^ Y}n{u^l + e|) (2.33) 



where e| = fc^ + is the energy of a one particle state with a certain momentum. 
Summation in respect with a;„ in ( |2.33| ) may be done with the help of the basic 
equation 



n=-oo ^ "T 'i ^ 
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The derivative of the sum gives 



2 exp Pe - 1 



(2.35) 



After integration with respect to e we will have 



2(3 



| + -^ln(l-exp[-/3e]) 



(2.36) 



The expression for the thermodynamic potential will have the standard form 



(2.37) 



Here the infinite, temperature independent part of the function Zj^ is dropped 



2.3 Green's function of boson field 



Introduce the two point Green function (or propagator) of boson field at finite temper- 
ature as the thermal average of a time-ordered product of two scalar field operator^. 



^In quantum field theory two-points Green's function is introduced as the vacuum expectation 
value of a time-ordered product of two field operators. 
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= Tr [exp{-pH)Tip{x)cp{y)] /Trexp(-/3i/) 



(2.38) 



One may express the Green's function as the sum 



Gp{x - x) = 9{t - T')Gt{x - x) + 9{t' - t)G-Jx - x) (2.39) 



where are Wightman functions 



G^{x — x!) —< Tip{x)(p{x') >p 



Gp{x — x) —< T(p{x')(p{x) >p 



(2.40) 



and 



eir) = 



1,T>0, 

[0,r<0. 



is the step function. 
In the interval [0, P] 



and 



(2.41) 
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(2.42) 



Using the fact that exp{—(5H) and time ordering operation commute and taking into 
account the cychc property of the trace, one gets an important property of the thermal 
Green's function 



Gpix - y)\x*=o = G^(x - y) 1^4=0 
= IV [exp{-pH)if{x, rMy, 0)] /Tr exp(-/3i/) 
= Tr [ip{y, 0) exp{-pHMx, r)] /Tr exp(-/3i/) 
= Tr [eM-PH){exp{PHMy, 0) exp{-(3H)Mx, r)] /Tr exp(-/3i/) 
= IV [eM-PHMy, PMx, r)] /IV exp(-/3i/) 



This relation leads to the Green's function in Euchdean quantum field theory. 



^Gp{x- y)\xi=i3 = Gi3{x - y)\xi=/3 



(2.43) 



Thus one has a periodicity condition 



Gpix - y) 1^4=0 ^G/six- y)\x4=p 



(2.44) 
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Since the Green's function Gp is periodic ( |2.44| ) and depends only on coordinate 



differences Gp{x — y ,x — y), it may be represented by Fourier series and integral as 



Gfi{x-y) = (l//3)^exp[z^„(x^ - y^)] J -^-^ exp[iA;(f 



(2vr)= 

xG^sK,^) (2.45) 



with ujn = 27m//3. 

In compact form one can write this equation as 



Gp{x -y)= eM^k{x - y)]Gp{k) (2.46) 



k 



where 



(2.47) 



The Green's function satisfies the equation 



(-dl + m^) Gp{x ~y) = 5{x - y) (2.48) 



with periodic conditions ( |2.44| ) and 5 -function in the form: 



—-^exp[i{u;nx' + kx)] (2.49) 



(27r) 
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From ( p:i6D and (|^) it follows that Gpik) is 



^2 _|_ jjj2 



(2.50) 



and 



exp{ikx) 

f^2 _|_ 



(2.51) 



2.4 Notation 

To work with the standard quantum field theory of signature (+2), return to imag- 
inary time formalism, on the time variable t = —ir defined in the interval [0, —ijS]. 
Therefore the equation ( |2.44| ) may be rewritten as 



Gf^ix - y)\xO=o = Gf3{x - y)\xO=-if3 (2.52) 
It leads to the replacement a;„ iui^ and 

//W = ^i:/(gl/K.*) (2.53) 

k 

with ujn = 2nin/(3. 

The equation for the Green's function becomes 
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k 

This expression looks like the expression for the Green's function in a common field 
theory. 
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2.5 Partition function for fermions 

The previous section discussed a quantization scheme for Bose fields in the functional 
integral formulation. However, the methods of path integral may be applied in the 
same way to finite temperature Fermi systems. This section develops these methods. 
Consider a free spinor field which is described by the action 

S = - J d'^xip{x){i-f ■ d + m)ij{x) (2.55) 

where {7/^} is a set of Dirac matrices, which are defined by 



{j„j.} = 2g,, (2.56) 

metric g^i, = diag{—, +, +, +)and, by notation, ip = ip~^j^ is hermition conjugate. 

The action (|2.55|) has a global U{1) symmetry and that is associated with con- 
served current 



f(x) = ip{x)Y'ip{x) (2.57) 



The total conserved charge is 



Q = J Sxf = J d^xtjj+i: (2.58) 



To apply the formalism developed in chapter II. 2 treat the field i/) as the basic field. 
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and the momentum congugate to this field 



dL 

d{dilj/dt) 



(2.59) 



Thus consider i/j and to be independent entities in tlie Hamihonian formulation. 
The Hamiltonian density is found by the standard procedure as 



dip I / dipX 

H = — L = tjj I ~ ^ ~ ipi—i'yo + m)ip 



(2.60) 



Introduce the partition function as 



Z = Tiexp[-ip{H - fxQ)] 



(2.61) 



Follow the steps leading up to the equation ( |2.21|) . For the partition function the 
intermediate equation will be: 



Za oc lim / idib^ — x 

^ N^oo J \ ^ 2tX j 



N 



exp <j -zAt ^ j (fx [H{%1)^, - ^fiipj+i - ^j)At 



(2.62) 



Taking into account the limit of this equation, find 



^f3 = N' J idt/j^ J dip 
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exphjdtj (fx li^jj^^ - H{i/j+{x,t),ilj{x,t))\ [ (2.63) 



For the time variable r, get 



exp |/ '^^ / ^^^'^'^ + ^7°7 ■ ^ - "^7° j (2-64) 

The quantization of a Fermi system can be obtained as a resuh of integration over the 
space of anticommuting functions iIj{x, r) (x G V^, r G [0, /?]), which are the elements 
of an infinite Grassman algebra. 

To obtain the correct result it is neccessary to impose on ip, ip the antiperiodicity 
conditions in the variable r: 

^(x, (3) = -^{x, 0), ^(f, P) = -ip{x, 0) (2.65) 

The Fourier series for i/j, ijj in Fermi case ( p.65| ) are 



V'(x,r) = {(3V)-'/'Y.T.^M^[^r + kx])Mk) 
" k 



ij+{x, r) = ipV)-^/^ Eexp(-^kT + kx])iP+{k) (2.66) 

" k 

—* — * 

where ipn{k) and tp^ik) are the generators of Grassmann algebra. 
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Inserting ( p.66| ) into ( |2.64| ) get 



Zh = N 



X 



J2J2'^n [{-i^n + /i) - 7°7 ■ - m7°] t/'^ (2. 



Calculation of the Gaussian functional integral over the fermi fields in (|2.67| ) leads 
the following equation 



Zp = Det[{-iuOn + /i) - 7°7 ■ k - m7°] 



(2. 



or using useful relation 



hiDetD = TrlnD 



(2. 



find 



In Zp = Tr \n[{—iuJn + fx) — 7°7 ■ k — mj^ 



2EEln [{uJn + ifif + e 



(2. 



Since both positive and negative frequencies are summed over, ( p.70 ) can be writ 
in the form 
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InZ^ = E E {in [(^n - f^f + 4] + In [(^- + f')' + 4] } (2-71) 



Summation (|2.71| ) over Un leads to the folowing equation 



In Zr- 



2^ E + (1 + exp{-/3(e,^ - /i)}) + (1 + exp{-/?(e,^ + /x)})] (2.72) 



The contributions for particles (/i) and antiparticles (— /i), and, also, the zero-point 
energy of vacuum, have now been obtained. [Kapusta 1989] 

2.6 Green's function of fermi field 

The finite temperature fermionic Green's function (or fermionic propagator) may be 
introduced in the similar way to scalar field. 
Consider that 



Ti[e-x.Y>{-(3H)T'tP{x),i){y)]/ Tr exp(-/3if) 



(2.73) 



For fermions the analog of ( p. 39 ) is in the form 



G^Jx -y) = e{x - y)G^^{x - y) + 9{y - x)G^Jx - y) 
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(2.74) 



where 



and 



Gg^ix -y) = -< iPiymx) < y' (2.75) 



The desired boundary conditions now foUow from 



Gff^ix-y) = -Tr[exp(-/3i^)V(|/^^)^(0,f)]/Tr exp(-/5i/) 



-Tr[exp{-pH) exp(/3//)V'(0, x) exp(-/3//)V'(?/^ y)]/TT ex.p{-pH) 



-Tr[eM-PHmP, ^i^{y\ y)]/T^ exp(-/3i/) = -G^^{x - y) (2.76) 



It leads to antiperiodic boundary conditions 



Gp{x - y)\x^=Q = -Gp{x - y)\xi=f3 (2.77) 



and this means that the fields ip, ip are antiperiodic 



iP{x,0) = -ip{x,l3) (2.78) 



The Fourier series of the fermionic propagator is written as 
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f 

Gp{x -y) = {l/P)^ex.Y> [iuJn{.x* - y"^)] J ^exp[ik{x -y)]G^{u;n,k) 

= Jexp[tk{x-y)]G^{k) (2.79) 

k 

with fc^ = {un, k), Un = (27r//5)(n +1/2). 

The Green's function of the fermionic field satisfies the equation 

{i^-d + m)G^{x-y)^S{x-y) (2.80) 



and has the following form 



Gj(^) = /^2^^-P(^M (2.81) 

k 



where {7;^} is the set of Euclidean gamma matrices. 

2.7 Notation 

As in the case of the scalar field one can work with a quantum model in Minkowski 
space-time, restoring the time variable t — —it defined in the interval [0, — It 
gives the following antiperiodic conditions for Green's function 

- y)|a;0=o = -Gp {x - y)\xO=-if3 (2.82) 
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and the rule of integration 




(2.83) 



with ujn — {2TTi/P){n + 1/2). In the Minkowski metric (signature (+2)) the equation 
for Green's function is written as [Dolan, Jackiw 1974] 

Gj(^)=/^^^exp(zM (2.84) 

k 
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Chapter 3 



THERMODYNAMICS OF 



QUANTUM GASES 



AND GREEN'S FUNCTIONS 

This section develops a method for calculation of thermodynamic potentials directly 
from finite temperature Green's functions. For this purpose the Schwinger proper 



time formalism [^chwinger 195 If has been applied, to write down generating func- 
tionals of quantum fields in terms of Green's functions of these fields. The finite 
temperature generalization of this formalism leads to thermodynamical potentials 
which are written through finite temperature Green's functions. 
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3.1 Thermal bosonic fields 



In standard form the generating functional Z of the scalar field >^ with the Lagrangian 
(|2.26|) is written as : 



Z\J\ = J L)(^exp|-(z/2) J d^xip{x)K^yip{y) + i J d^xJ(x)(^(x)| (3.1^ 



The symmetric operator 



K,y={-dl + m')6{x-y) (3.2) 



can formally be treated as a symmetric matrix with continuous indices {x,y). It has 
the following properties: 



X — z) 



I d'vKi'Kli' = K^, (3.3) 

The functional Z gives the transition amplitude from the initial |0~ > and final lO"*" > 
vacuum in the presence of the source J(x), {Z =< 0+|0^ >j) which is producing 
particles. The Green's function may be treated as the solution of the equation 



(PK^yGiy, x) = S{x - x') (3.4) 
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with operator K^y- Now rewrite ( |3.1| ) in a convenient form changing the integration 
variable from ip to 



ip\x) = / d'^yK^yip{y) 



(3.5) 



Then the quadratic form of (3J.) may be recast as 



{1/2) jd'x \v\x)- j dSK-^'i'J{y) 



{1/2) dVy^{y)G{x,y)^{y) 



(3.6) 



Substituting ( p.6|) into (|3.1| ) and performing an integration of the Gaussian type 
integral, get 



Z[J] oc (deti^^/2)-^exp 



i/2) d''xdy^{y)G{x,y)^{y) 



(3.7) 



where the Jacobian arises from the change of variable ( ^.51 ). 

The functional determinant in ( |3.7| ) may be written with Green's function G(x, y): 



{DetK^'^)-^ = {DetG{x,y)y/'' = exp [(l/2)TrlnG'(x 



1/2 



(3.8) 



The Green's function is found from functional differentiation of Z with respect to 
source J 
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Now pay attention to the functional determinant in ( p.7|) because of its important 
role in the applications to statistical mechanics. 

One can introduce the heat kernel for operator K as the solution of the partial 
differential equation 



d 

i—^{x,x;is)= I d'^zKxz^{z, X -jis) (3.10) 



with the boundary conditions x'; 0) = 5{x — x) 

The Green's function G may be written with in the form 



oo 



G{x^x) = J ids^{z, X ;is) (3-11) 



and the logarithm of functional determinant of the operator K as 

oo 

lnDetK = J ids{is)~hTQ{z,x ^is) (3.12) 



To get the equation for the heat kernel 3, find proper time representation of Green's 
function. 

Using the useful equation 

oo 

{k^ + m^)-^ = i J dsexp{-isik'^ + m^)} (3.13) 
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write the equation for Green's function in the form 



(27r)4 P + m? 



J ds J exp{iky — is{k'^ + rn^)} (3-14) 







(2vr) 



where y = x — x . 

After the integration with respect to momentum, we get the Schwinger represen- 
tation of the Green's function 



oo 



G{x,x) = -——^ J idsiisY eyi-p^—isrn? — a /2is) (3.15) 







{Any 

The variable a equals half the square of the distance between x and x 

a = {x-x'f/2 



The equation for the heat kernel directly follows from the comparision (|3.11 ) and 



^{x, X \ is) = -r exp(— ism^ — a I2is) (3.16) 



For future calculations it is better to work with the generating functional of connected 
Green's functions W[J] which is connected with generating functional Z by means of 
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equation 



Z[J] = exp{ziy[J]} 



(3.17) 



From ( p:T2D and (|3T7|) we find 



W[Q] = {i / 2) In DetK 



{i/2) J ids{is) ^tr53(x, x ; is) 





(3.18) 



Inserting ( p.l6| ) into ( p.lSj) and taking into account ( p.l5|) we get an important equa- 



tion 



oo 

4„ / J 2, 



W[Q] = -{i/2) / d'^x / dm'tTG{x,x) 



(3.19) 



Now find finite temperature functional inserting Green's function at finite tem- 
perature 



oo 

j4_ / J 2. 



= -{i/2) J d X J dm^iiGi3{x,x 

/3 m? 



(3.20) 



Now 



Wp = -(/3/2) J d^x J dm^ j {e + m^)-^ 

V rrfi k 
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= i/3/2)V J\n{k^ + m^) (3.21) 

k 

The Helmholtz free energy may be treated as effective potential of tlie finite part of 
Wjs, therefore 

V) = -iWp. (3.22) 



and 



F{P, V) = -{iP/2)V J + m^^ 



= (l/2)l^/^EM^n + 4) (3.23) 

After making the summation over the frequencies u, this equation will coincide with 
( |2.37| ) (if the infinite, temperature independent term is dropped). 

So one can conclude that the density of Hemholtz free energy may be written as 
[Kuhkov & Pronin 1987] 



f{f3) = {i/2) f dmhTG{(3,x-x'), (3.24) 



where G{(3,x — x') is temperature contribution in Gf3{x,x] 



3.2 Bosonic finite temperature Green's 

function in the Schwinger representation 
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Start from the finite temperature Green's function ( |2.51| ) in the hmit of coincidence 



X = X : 



hm^ Gplx, X 



d^k 1 
(27r)4 fc2 + m2 
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/ ds{l/P) E / exp{-^s(-c.^ + P + m^)} (3.25) 



After the integration with respect to momentum find 



00 

\im^Gf3{x,x') = /"(47ris)3/2gxp(_isr7i2)(i//5)^exp(-47rV//32), (3.26) 



Rewrite the sum in ( |3.26| ) with the help of the equation ||Ditrich 1978 



Eexp[-a(n -2)2] = Y.i'^/aY^^ exp{-7i^n^ /a - 2nizn) (3.27) 



In bosonic case {z = 0) ( p.26|) has the form: 



47r r „ J us r 



(3.28) 



Selecting the temperature independent (n = 0) part of the equation (|3.28|) we find 
the finite temperature contribution in Green's function of scalar field: 



Gi3{x, x) = G{x, x) + G{P, X — x) 
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i r ids 



(47r)2y {is 



% ^ — ^ f %ds . 9 9 /n9 / \ 

(4vr)2 (^s)2 



— expi^— «sm 



(3.29) 



Inserting G{[3,x — x) of the equation (|3.29|) into ( 3.24 ) we find the equation for free 
energy in the form 



7 °^ 7 ids 
/(/?) = -(l/(47r)2) dm^Y. 7^exp(-2sm2 -n2/?V4«s 



(3.30) 



The integral over (s) occuring in (|3.3(JD can be found with 



j dxx"'^ exp (^-a- - 7^^ = 2 (a/7)''/^ ir^(27^), (3.31) 



where is a modified Bessel function. Then 



ids{isy~^exp{-ism'^ - P^n'^/Ais) = 2 (/?n/2m)-''~^ Kj_2iPmn), (3.32) 



and the equation for / will be 



n=l 



(3.33) 



The integral representation of the sum of the modified Bessel function ( PCXII.32"D 
allows the standard equation for the Helmholtz free energy to be written as 
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/(/5) = (l//5)/ (^ln(l-exp(-/?e,)) 



(3.34) 



directly from (|3.33|) . 



3.3 Thermal fermionic fields 

Let us now develop the same formalism for fermionic fields. 

Add to the fermionic Lagrangian the terms with anticommutative sources r]{x) 
and f]{x) 



L — > L = Lq + ri{x)i/j{x) + 'ip{x)ri{x) 



(3.35) 



and write the generating functional 



Z[ri,T]]= / DipDijj exp li / d'^x Lq + ■r]{x)'ilj{x) + ip{x)ri{x) >, (3.36) 



where the functions ip and ip are considered as Grassman variables. 
After integration in respect with these Grassmann variables obtain 



Z[ri, rj] Qc {DetK) exp 



-i J d^x J d'^yri{x)K^yr]{y) 



(3.37) 



where 



Kxy = {il ■ + m)5{x - y) 



(3.38) 
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and bi-spinor Gp satisfies the equation which is equivalent to 
In operator form [[Schwinger 195 1|| it is written 



KGp = 1 



At zero sources the generating functional W[0] is 



W[0] = ~i\nZ[0] = -t\nDet{K) 



On the other hand it is 



W[0] = -i\nDet{K) 



{i/2) J idsiis) ^ti^{x,x]is) 





where the kernel '^{x, x' ; is) is the solution of the equation 



d ^ , - / 

i—Q{x, X ^is) = {KQ){x, X ]is), 



The equation for the bi-spinor Green's function is written as 



Gf{x, X 



d^k 1 
(27r)4 P + m- 



exp{iky) 



67 



ds J j^-^exp{iky - isik"^ + m'^)} 



(3.43) 



and after integration over the momentum, we get 



oo 

/ 1 r 

Gf{x, X ) = — / dsiisY exp{—ism'^ — a/2is) 

(47r)^ J 



(3.44) 



As follows from (|3.11|) the heart kernel has the form 



Q{x, X ]is) = — Y ^^^)"^ exp{—ism'^ — a/2is), 



(3.45) 



Inserting ( |3.45|) into ( p.41| ) and using ( |3.44| ) we find 



oo 

W[0] = (z/2) / d'^x [ dmhTGF{x,x) 



(3.46) 



For non-interacting fermionic field this equation is divergent and does not give any 
useful information but, in the finite temperature case we can get some interesting 
results connected with temperature properties of quantum Fermi gas. 



3.4 Fermionic finite temperature Green's 
function in Schwinger representation 

In the limit of coincidence {x x') the finite temperature fermionic bi-spinor Gp 
is written as: 
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Im Gf{x, X ) = 



d^k 1 

(27r)4 P + m2 



/ E / exp{-z.(-^^ + k' + m')} (3.47) 



(27r) 

where u;„ = (2i7m//3)(n + 1/2). 

After the integration in respect with momentum we get 

hm Gpix, X ) = 



/ {4msf/^exp{-tsm^){i/(3)J2exp[-A7r\n + l/2)yf3^] (3.48) 



Using the equation ( p. 271) with z = —1/2, we find 



^ exp[-(27r//5)2(n + 1/2)^] = ^(-l)"[/5/(47rzs)^/2] exp(-7r2nV4zs), (3.49) 

n n 

Then the equation ( |3.48| ) will be 

lim Gix, x) = i-J-^ y (-1)" J ^exp(-ism^ - n^p'^/Ais). (3.50) 
Now we can select the temperature independent part of the equation (|3.50|) : 

Gp{x, x) = Gf{x, x) + Gf{P, X — x) 
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i r ids 



(47r)2y {is 



— expi^— 2sm 



+ ^ £ (-1)" / 1^ eM-^sm^ - n^PV^^s) 



(3.51) 



Inserting (|3.51|) into the equation [Kulikov & Pronin 1987]: 



f = {-i/2) / dmhiGpifi^x- x)) 



(3.52) 



we find an expression for tlie free energy of Fermi gas in tlie form 



OO rl 

f(l3) = -(l/(47r)')tr(i) / dm^ J2 ("l)" / 7-^2 expi-ism^ - n'p^/Ais) (3.53) 



The integral over (s) occuring in ( |3.53|) can be written as 



ids{isy-^ exp(-ism2 - /?V/4is) = 2 (/?n/2m)^"^ Kj_2{pmn), (3.54) 



and the equation for density of free energy of fermionic field / will be 



00 ( 1 An 



/(/3) = tr(mV2vr2/32)(i) ^ ^ir,(/3mn 



(3.55) 



The standard equation for Helmholtz free energy follows from ( |XX11.1^ 



/d h 



(2vr) 



(3.56) 
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Now we can describe the thermal behavior of massless vector fields in the Schwinger 
proper time formalism to prepare the mathematical foundation for further computa- 
tions in curved space-time. 
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Chapter 4 

FINITE TEMPERATURE 
GAUGE FIELDS 



In the present day progress in quantum field theory is to a great extent due to the 
development of gauge fields [[Yang fc Mills 195^ , [Utiyama 1956| , ft Hooft 1971|| . These 



fields open up new possibilities for the description of interactions of the elementary 
particles. Gauge fields are involved in most modern models of strong, weak and 
electromagnetic interactions. They are extremely attractive for the unification of 
all interactions into a single universal interaction. On the other hand the functional 
formulation of the gauge fields helps us to get statistical and thermodynamical results 
connected with finite temperature properties of the particles which are described by 
these fields. In this connection we will develop the finite temperature formalism for 
gauge fields 
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4.1 Gauge theories: Pure Yang-Mills theory 

The basic idea of gauge field relies on the local gauge invariance principle of the 
quantum field theory |[Y"ang fc Mills 1954| , [Utiyama 1956|| . Let us consider a gauge 



tansformation parametrized by the functions (^"-{x) 



U{x) = exp[-icj"(x)r'^] (4.1) 

where the generators of the Lee algebra obey the equation 

[r", = ir^'r" (4.2) 

Numbers /"^^ are the structure constants of the group. 

Let the vector (gauge) fields Afj^{x) = A^(x)"r" transform according to 

A^{x) A;{x) = UA,{x)U-' + {i/g)Ud,U-^ (4.3) 

For infinitesimal transformations we have 



{A^l = + r"'Al{x)uj\x) + {l/g)d,u;\x) 



Al{x) + 5Al{x) (4.4) 
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The transformations (|0 



) preserve the Maxwell type Lagrangian 



L = -(l/2)Tr(F,.(x))2 = -(1/4)F;,F'^'^'^ (4.5) 

where 

Ff,^ = df,A^- d^A^ - ig[A^, AJ\ 

= {d.At - d^Al + gr'''A\Al)T'^ (4.6) 

is the strength of the vector field A^. 

The path integral over the field A^ with the Lagrangian ( |4.5| ) in the form 

1 DA^exp(z5[A^,9,A^]) (4.7) 

is undefined because of the gauge degrees of freedom. Namely, the gauge transforma- 
tions create the "orbit" of the field A^^ in the functional field space and the functional 
integration (|4.7| ) over this space overcounts the degrees of freedom of the theory. To 
improve the situation we have to eliminate all unphysical degrees of freedom which 
origin from local gauge invariance. So, we have to "slice" the orbit once so that we 
do not have this infinite overcounting. This is the origin of the gauge fixing problem. 
To solve this problem one can select the surface in this functional space. The surface 
is good if it intersects the orbit of any given field under gauge transformation once 
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and only once. The equation of this surface may be written as 



n^,) = 



(4i 



To ehminate the unphysical degrees of freedom we may insert the factor ^[^(A^)] in 
the functional integral. 
We can do it inserting 



1 = A^p(A^) j Du5[F{A^^)] 



(4.9) 



into ([4.7|) . This expression can not change the measure of integration in ( [4. 71) so we 
will have 



j DAfj_ exp |i j d'^xL 



J DA^ (App(A^) j Duj6[F{A;)]^ expj^ J d^xL 



(4.10) 



where 6fp is Faddeev-Popov determinant, and Duj = Y[Y[duj"'{x) is the invariant 

a X 

group measure fFaddeev fc Slavnov 1980|| . 

One can show that it is gauge invariant fRyder 1985 



Afp(^m) = ^Fp(^: 



(4.11) 
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therefore the equation ( [4 .101 ) is 

I DA^j D^A^p(A-)5[F(A-)]exp|^ j d'^xL^ 



(4.12) 



Making the gauge transformation from to A^ and taking into account that action 
is gauge invariant too, we get, dropping multiphcative divergent factor / duj, the 
expression for generating functional of the vector field A^ 

Z[J] = j DA^A^p(A^)5[F(A^)]exp|z j d^xL + i{J,A)^ (4.13) 

Our task now is to calculate an explicit result of the Faddeev-Popov determinant 
App(A,J. 

To do this we denote the numerical value of the function F{A^) at x as 



F{A,) = F\x) 



(4.14) 



Then one can write the equation (4J) in the form 



^-f'p{A) =\{\{ j du\x) 5{F\x)) 



Xm dF\x) 5{F\x)) 



d{uJi{x)...ujN{x)) 
d{F,{x)...Fj,{x)) 



n 



duj'^ix) 



dF\x) 



F=0 V"-^ / |F=0 



(4.15) 
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The Faddeev-Popov determinant 



AMA) = Det { ^] (4.16) 



is the functional determinant of the continuous matrix || 5F"-[x) / 5uj^{y) || with rows 
labelled by (a, x) and columns by (6, y). 

The generating functional ( [4 .131) will have the form 



Z[J] = J DA^{x)Det (^SFySu^) S[F%A^)] exp {iS + i{J, A)} (4.17) 

Now we can make the last step to rewrite our expression in a form convenient for 
practical calculations introducing Faddeev-Popov (ghost) fields. 

4.2 Ghost fields 

For the calculation of generating functional ( [4.17[ ) we must select a gauge constraint 



and compute the FP-determinant ( [4.16|) in this gauge. 



Let us select the Lorentz gauge constrain 

F'^(A^) = d^'A^^ix) = (4.18) 

We can find Det{5F^/5uj^) from the Teylor expansion of the with gauge transfor- 
mation uj^: 
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pa _^ ^pu>Y = pa^ Det{5F''/5uj^)uj^ + ... 



(4.19) 



For our gauge conditions ( [4.18|) it is 



d^Al{x)^d^Al{x) + r''^d{Al{x)u:\x)) + (lA)9V(x) 



d^Al{x) + {{l/g)dH'''^ + Al\x)d^)uj'^{x) 



d^Al{x) + jdyY. [{l/g)d'5'''5{x - y) + Al\x)d^5{x - y)\ u\y) (4.20) 



where A"^ = /'^*'=A°(x). 



Then we can find, that 



Det [SF/6uj] = Det \\< x,a\SF/Suj\y,c>\ 



Det 



{{l/g)dH'^^ + A^^{x)d>^)5{x-y) 



(4.21) 



One can rewrite Faddeev-Popov determinant (|4.21| ) in terms of a path integral over 

anti-commuting fields (c, c)[]. They are called Faddeev-Popov ghosts, the fields c, c 

are independent anti-commuting scalar fields which obey Fermi statistics. 
^The functional determinant of anti-commuting fields c, c is 



/ 



DcDceK^{i^M''^J'} = Det{M) 
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Then the determinant ( ^.21| ) may be written in the form 



Det[6F / 6uj\ = J DcDcexp i J d^xLghost{x) 



(4.22) 



Lagrangian Lghost{x) is 



■'ghost 



X 



c''{x){d^S'''' + gAl\x)d'')c\x) 



(4.23) 



The first term in the Lagrangian ( ^.23| ) is the pure ghost part 



(4.24) 



and the ghost propagator Go directly follows from ( 4.24 ). The equation for this 



propagator is 



dlGo{x-y) =S{x-y), 



(4.25) 



In the momentum space Go{k) is written 



The second term in (|4.23| ) describes interaction ghosts and the vector field 



(4.26) 



For better understanding of the ghost contributions let us study an effective action 
of the ghost fields. 
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4.3 Effective action 



The functional determinant (|4.21|) may be rewritten as 



Det{6F/6uj) = Detid'd"' + gAl^d, 



Det{{l/g)d^)Det(l + g{l/d')Al'd, 



Ignoring the first determinant^ 



Det{d^) = Det{GQ^), 



we can write an effective action in the form 



Seff = S -iTrln{l + L), 



where the element L is written as 



< X, a\L\y, c>=g<x,b\ {l/d'^)Al^{x)di^Six - y) 



g I dzGo{x-z)A';^'ix)d^6iz-y). 



^This term is important in the finite temperature limit. 
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Now we can expand the determinant with 



Det{l + L)= exp {Trln(l + L)} 



exp y ^ '- TrL" (4.31) 



Write the contributions Tr(L") in ( ^.31| ). 
The first expression of the order g is 



\nL = tTjdxJ dz [Go{x ~ z)gA''^'{z)dJ{z - x)} (4.32) 



because tr(A;^^(x)) = tr(/'^^^A|^(x)) = 0. 

The second expression of the order will be 



In(L^) = (yf^tr J dx J dy J dz J dtGo{x — z)x 



X A'';{z)dJ{z - x)Go{y - t)Al%t)dt6{t - x) 



dz / dt {gA';^{t)dtGo{t -z)}x 



X {gA^;{z)d,Goit - z)} (4.33) 



As we can see from ( [4.33|) and ( [4.31| ), for > 2 the trace of ( [4.29|) equals the integral 
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with integrand in the form of the cycle multiphcation of the terms 



{^A^^(z)a,Go(t-^)} (4.34) 



namely, we can write 



Tr(L") = J dz,...J dzr,{gA';''{z,)d,,Go{z,-Z2)} 



X ... 



X ... {gA';%Zn)d,„Go{z^^i - (4.35) 

These terms give non-trivial corrections due to ghosts to the effective action S'ej/[yl] 
Therefore the determinant (|4.31| ) may be treated as contributions of the closed 



loops with internal ghost lines, and external lines of which correspond to vector fields. 



4.4 Propagator of vector field 

In the previous section we found out how to compute the Faddeev-Popov determinant 
in the generating functional (|4.17| ). Now we must treat the second important term 



the delta function arising from the gauge constraint. For this we will rewrite the delta 
function as [[Vladimir ov 1976[] 



sid'A';,) = Y[6id'^A';,ix)) 



X 
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lim 2i7ra) "'^''^exp 



2a 



~ lim exp < I d X 



(4.36) 



The argument of the exponent ( ^.36| ) may be combined with the quadratic part of 
the vector Lagrangian 



i d X 



(d'Alix) 



(4.37) 



The second term of the expression ( [4. 3 7| ) fixes the gauge. 

As a resuh we have the quadratic part over the vector fields in the form 



\jd'xAl{x) 



Atix) 



(4.38) 



The propagator of the vector field may be found from the equation 



G,xix~y) = -S,xSix-y). (4.39) 



The solution of this equation is 



d^k 



rjf,^ + (1 - a)-p- 



exp(iA;x) 
k'^ + ie 



(4.40) 



The finite temperature propagator of the vector field in the Feynman gauge (a = 1) 
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has the simple form 



GM= / j^^eMikx) (4.41) 



k 



where the definition of the integral / is the same as in ( |2.54|) 

k 



4.5 Partition function for Gauge fields 

Let us study thermal gauge fields in the axial gauge (A3 = 0). 
The standard expression for the partition function is 

Z = Trexp[-PH] 

where H = Hlnf, A^) is the Hamiltonian of the Yang-Mills field. 
According to section II (|2.15|) we can write 

Z = i^Texp[-pH])a,^al (4.42) 



periodic 



X exp y c/r y d^'xlm.AI^ - A^)] | (4.43) 

The Hamiltonian in this gauge has only two degrees of freedom for each vector field: 
^4^2- The conjugate momenta for these fields are tii2- 
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After Gaussian integration over the momenta, we get 



Z = N 



DA'' X{5[Al] exp Ijdrj d^xL{x) 



periodic 



(4.44) 



The result is similar to the expression for the Faddeev-Popov ansatz in the axial 
gauge. Calculations of ([4.44|) give the correct result for the partition function of a 
photon gas. The same result will also be obtained for Coulomb gauge (VA = 0). 

For the Feynman gauge the simple calculation of the type ( [4.43| ) leads to an 
incorrect result. But we may correct the result taking into consideration the ghost 
contribution and rewrite the partition function ( [4. 431 ) in the form of ( |4.13| ). 



Z = N J DA''Det{5F''/5uj^)5[F''] exp \ J dr J d^xL{x) 



periodic 



(4.45) 



Since under the gauge transformations = —duU"", we have 



Det{6{df,A^'')/6u^) = Det{d^) 



(4.46) 



that coincides with ( |2.28| ) at order ~ ^{g)- For the Feynman gauge (a = 1) 



DA" exp \ / rfr(-l/2) / d^x \Al{x)d^ A"^^''^] \ = Det{G^,y/^ (4.47) 



periodic 



At finite temperature Det{G^u) may be written in the form 



Det(G'^,) =exp(-l/2)Trln(G'^,) 
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= exp(-2)X:/^ln(a;^ + ^^) (4.48) 



and the temperature contribution of Det{d'^) in the form 



Det{d'^) = expTrln(a^) = expTrln(G) 



expE/(0jM^^ + ^') (4.49) 



where G is just the propagator of the scalar field. 

The logarithm of the finite temperature generating functional 2^[0], according to 
previous results, is 

lnZ[0] = \nDet{G) - {l/2)\nDet{G^^), (4.50) 



or 



ln^[0] = -E/^ln(a;^ + ^')- (4-51) 



After summation the free energy is found in the form 

/(/?) = -^lnZ,= 0) / _g_ln(l-expM6]), (4.52) 
where = {ky. 
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So, we got the correct answer || Bernard 1974|| by introducing the finite tempera- 



ture ghost contribution. The index (2) in the equation ( 4.52|) refiects the two degrees 



of freedom of the massless vector field. The generahzation to Schwinger proper time 
formahsm is the same as in the previous section. 



87 



Chapter 5 



QUANTUM FIELDS 



IN CURVED SPACE-TIME 



In this chapter we introduce the basic formahsm of quantum fields in curved space- 
time which was studied in works of DeWitt, Fulhng, Bunch et al. fPeWitt 1975|| , 
Fulhng 1989^ , [Parker 1979|| , jBirrell fc Davies 1982 1. 



5.1 Lorentz group and quantum fields 



Field theory in Minkowski space-time has been well studied from the group point 
of view by many authors [ Sterman 1993 |,| Ramond 1983 |. Here we will consider the 
connection of fields of different spin with the Lorentz group. 

In fiat space the spin of the field is classified according to the field's properties 



under infinitesimal Lorentz transformations 
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= Apf^ = (5^ + uj^)x'^ (5.1) 



with 



<^a/3 — <^/3a 

which preserve the length of the coordinate vector x'^. 

Under Lorentz transformations the general multicomponent field 

pap...\ 



transforms according to 



^a/3...A ^ [D{A)fJ^;\,F''^ (5.2) 



where 



[D{A)] = 1 + {l/2)u;^^E^p. (5.3) 

In order for the Lorentz transformations to form a group, the antisymmetric S^/j 
group generators are constrained to satisfy 



pa/3, S^i] = rj^p^aS — Va-Y^pS + VS/S^'ya — VSa^-yP 
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(5.4) 



We may write down Sq,^ for different types of fields: 

1) It is easy to see that the scalar field has the following form of S 



= (5.5) 

2) A vector field transforms as 

^ A;^,F"' (5.6) 
so from (|5.3| ) and ( p.l| ) we get an expression for S in vector case in the form 

\^c,pV& = ^IVf^S - ^JVaS- (5.7) 

3) For a spinor field it may be written as 

[S„^]^, = (l/4)[7„,7/.]^„ (5.8) 

where 7 are Dirac matrices. 

Now we can generalize this concept to curved space-time. 

5.2 Fields in curved space-time 

Let us consider general coordinate transformations. We will define a general coordi- 
nate transformation as an arbitrary reparametrization of the coordinate system: 
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= x'^x") (5.9) 



Unlike Lorentz transformations, which are global space-time transformations, general 
coordinate transformations are local. Under reparametrizations, a scalar field trans- 
forms simply as follows: 

^'{x)=^{x), (5.10) 



and vectors and dx^ as: 



dx'^ 

dx^ = — — dx'^ (5.11) 
dx" ^ ' 



As we can see, the transformation properties of (|5.1CI| ) and (|5.11| ) are based on the 



set of arbitrary real (4 x 4) matrices of the group GL(4). Now we can give the 
abstract definition of covariant and contravariant tensors and write the general form 
of transformation 



= |^...|^...4V::(-) (5.12) 



Let us introduce a metric tensor g^^^ which allows us to write the four interval 



ds^ = g^ydx^dx 
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(5.13) 



The tensor g^j^i, transforms under the general coordinate transformation as genuine 
tensor. The derivative of the scalar field d^iplx) is a genuine tensor, but the derivative 
of a vector is not. To create a vector from dij,A,^{x) and {dij,A'^{x)) we introduce new- 
fields, connections, that absorb unwanted terms. The covariant derivatives are written 
then as 



V^A^ = d,A, - r;,A^ (5.14) 

Prom the restriction 



Vg^, = (5.15) 

we get the expression for connection Fj^^: 

rj, = {l/2)g^''{d^g,, + d,g^, - d,g^,) (5.16) 

which are named the Christoffel symbols. 

Another very important object in curved space-time is the Riemann curvature 

tensor R^^j,^^ which arises from the commutator of the covariant derivatives 

[V^,V.]>1^ = 4,,A^ (5.17) 
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The Riemann tensor is written in the form 



Contracting the indicies of the Riemann tensor one can get the Ricci curvature tensor 



R,. = RU (5.19) 



and scalar curvature 



R = g^'R^.u (5.20) 

The transformation properties of the volume element and the square root of the metric 
tensor are 



d^x = Det ( ^] d'x 



g'(x')=Det(^)^) (5.21) 



From ( ^.211) we get that the product of these two is invariant 



g{x)d^x = inv. (5.22) 
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Now we can construct the Einstein-Hilbert action as an invariant of the form 



Sg = -{l/2k') J d^xy/^)R, (5.23) 



and the action for a scalar field 



= -(1/2) J ^) (g^'-d.ipd.ip + mV) , (5.24) 



where the scalar matter couples to gravity via the interaction 



The coupling of the gravitational field to the vector field is also straightforward 



Sv = -(1/4) / V^(x)^''V«F;,F« (5.25) 



For the gauge-fixing term the generalization to curved space-time is expressed by 



S,f = -{l/2)a-' J d'x^) {V,An' (5.26) 



For the ghost term it is of the form: 



Sgh = / ci^xy^(-a^c"a'^c") (5.27) 
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However, the couphng of gravity to spinor fields leads to a difficulty, because there 
are no finite dimensional spinorial representations of GL{4). This prevents a naive 
incorporation of spinors into general relativity. The method we may use for these 
constructions involves vierbein formalism. 

5.3 Spinors in General relativity 

The vierbein (tetrad) formalism utilizes the fact that we can construct a fiat tangent 
space to the curved manifold and introduce the spinorial representation of the Lorentz 
group in each point of the manifold in this tangent space. Spinors can then be defined 
at any point on the curved manifold only if they transform within the fiat tangent 
space. 

Let us erect normal coordinates y^^^ at each point X of the space-time manifold 
M. To preserve the connection with the previous sections we will label the fiat tangent 
indices with letters a, /3, 7, 5,... from the beginning of the Greek alphabet, as we 
did it before, and general coordinate transformation indices with letters A, /i, u, 
from the end of the Greek alphabet. 

Introduce the vierbein as the mixed tensor (matrix) 



Note that the label a refers to the local inertial frame associated with normal coor- 




a = 0,1,2,3 



(5.28) 



dinates y^^^ at the point X, while // is associated with the general coordinate system 
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The inverse of this matrix is given by 



= S'^'^ (5.29) 



The vierbein can be viewed as the "square root" of the metric tensor gfj^^ : 

g,u = hy^v^p (5.30) 

For general coordinate transformations x'^ = x^{x^') we can consider the effect of 
changing the x^ while leaving the i/p^-) fixed. Then the verbein transforms as 



la OX 



hl-^h= (5.31) 



We also can transform the y'^-^^ arbitrarily at each point X 



yfx) - yjx) = Kx)Xx) (5.32) 

In this case h'^{X) transforms as a Lorentz covariant vector 

K(X) - ^1{X) = A{X)XiX) (5.33) 

which leaves the metric ( ^.30| ) invariant. 

If a covariant vector is contracted into h'^, the resulting object 
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A, 



a 




(5.34) 



transforms as a collection of four scalars under a general coordinate transformations, 
while under local Lorentz transformations ( |5.1| ) it behaves as a vector. 

Thus, by use of tetrads, one can convert general tensors into local, Lorentz- 
transforming tensors, shifting the additional space-time dependence into the tetrads. 

Now we may construct the generally covariant Dirac equation. 

We introduce a spinor il}{x) that is defined as a scalar under general coordinate 
transformations and an ordinary spinor under fiat tangent space Lorentz transforma- 
tion: 

Coordinate transformations: ip ^ ip 
Lorentz transformations: ijj — D[A{x)]ijj 

It is important to note that we have introduced local Lorentz transformations 
in fiat tangent space, so uJa/B is a function of the space-time. This means that the 
derivative of the spinor is no longer a genuine tensor. Therefore we must introduce 
a connection field u'^^ that allows us to gauge the Lorentz group. The covariant 
derivative for gauging the Lorentz group may be written as 



(5.35) 



Let {7"} be a set of Dirac matrices with 



{7a, 7/3} = 277„/3 



(5.36) 
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in the tangential space-time. 

The Dirac matrices 7" can be contracted with vierbeins: 



K{x)l. = l,{x) (5.37) 



Then 



{l^,{x),-i^{x)} = 2g^^{x) (5.38) 



In the result ( |5.33| ) and ( [5.35| ) the generally covariant Dirac equation is given by 



{i-f^{x)V^^, + m)i,{x) = Q (5.39) 
and hence the action for Dirac particle interacting with gravity is given by: 

L = ,/^ilj{x){iY{x)V^,,^ + m)ilj{x) (5.40) 

where ^/g = Det{h^). 

We can construct a new, alternative, version of the curvature tensor by taking the 
commutator of two covariant derivatives: 



[V^,V.]^ = ^i?^^S„^^ (5.41) 
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Written out, this curvature tensor is generally covariant in v, but flat in a, f3: 



= d,u;f - d^ul^ + ul^ul^ - uZ^u^f (5.42) 

At this point, the spinor connection o;";^ is still an independent fleld. Covariant 
derivative of the object with two different kinds of indicies is written as: 

^^Al = d,AZ + V^A^^ + uj\pAl (5.43) 

The external constraint 

= W + ^UK + = (5.44) 
helps to express the spin connection through veirbeins 

+ '^h'^^'h^^^d^K^ - d,h^,)hl -{a^/3) (5.45) 

After this preliminary work we can study properties of bosonic and fermionic fields 
in an external gravitational field. 
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Chapter 6 



BOSE FIELDS 

IN CURVED SPACE-TIME 

As we know already the inclusion of the interaction with gravitational field (theory 
formulated in curved space-time) is accomplished by replacing of the partial derivative 
by the covariant derivative 9^ — ^ V^. It is necessary to ensure that the Lagrangian 
is the scalar under the general coordinate transformation. Integration is performed 
over the invariant volume. This procedure, based on general -coordinate covariance, 
is called the minimal interaction for gravity and leads to the action ( ^.24 ). 

However, general coordinate covariance does not forbid adding to the Lagrangian 
invariant terms which are vanishing in fiat space-time. Such terms describe the non- 
minimal interaction with gravity. Therefore the theory under such consideration can 
be written in the form 
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S = J d^X\/ —g {L{(f), Vfj,(j)) + non — min. int.) (6.1) 

From dimensional analysis of R and conclude that the term, describing non-minimal 
interaction of matter field with gravitational field may be written in the form: 

(1/2)^V (6.2) 

where ^ is a dimensionless parameter (non- minimal coupling constant) ||Fulling 1989|| . 

If m = and ^ = 1/6 then the action is invariant not only under general-coordinate 
transformation but also under conformal transformations 
[Birrell fc Davies 1982 

V(^) = e'^^^^g,., ^'{x) = e^(^V(x) (6-3) 

where x(^) is an arbitrary scalar field (parameter of the conformal transformation). 
In the result we will have the action in the form 

= -(1/2) J d^x^^{x) + + iR) ^{x) (6.4) 

where the d'Alembertian operator is = giiu{,x)V^ fjS/^ u = d^d'^. 
The generating functional will be 

Z[R, J] oc {DetG)-^/^x 
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X exp|-(V2) J <fxg{xf'^j d^x g{xf'^ J{x)G{x,x')J{x')^ (6.5) 

where G{x,x) is the Green's function of the scalar field which is described by the 
equation: 

g^^'^\x) (-n^ + + iR) G{x, x) = 6{x - x') (6.6) 

So, for computation of the generating functional ( |6.5|) we have to define the Green's 
function of the scalar field in curved space-time from the equation ( |6.6| ). 

6.1 Momentum-space representation of 
the bosonic Green's function 

In curved space-time we cannot solve the equation ( |6.6| ) and compute the Green's 
function G(x,x') for arbitrary points x and x' of the manifold. But we can do so in 
the particular case of the limit of coincidence {x ~ x ). This limit gives us a possibility 
to find the (Det) of the Green's functions to get the effective action. Thus we will 
treat the problem of the Green's functions calculations in the limit of coincidence. 

Let us select a point of the space-time manifold x and construct the tangential 
space at this point as we did in the previous section. Let any point x of the manifold 
have the normal coordinate y°'{x). This coordinate may be treated as a vector in the 
tangent space with the origin at the point x . In this tangent space the norm of the 
vector y°'{x) will be (y,?/) = y'^{x)y^{x)riaf3-, and the metric of the manifold may be 
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written as [[Petrov 1969| 



(6.7) 



and 



9{x) = 1 - iRc^'-y^ - \Rap-ay''y^y'' 



^ ( T^^apR-yS - 77f.Ra(3xR^5u " 7^Raf3;-yS ) y"y^y'^y° + 



18 



90 



20 



cj /3 7 (5 



(6. 



where the coefficients are calculated at the {y = 0) origin of the coordinate system. 
The second derivative for the scalar field is 



y.^y. = v'^^d^dp + l-Rj^'y-y^dld^^ - ^i?^/^^ + ... 



(6.9) 



For our calculations it is convenient to express the Green's function G{x,x') as: 



G{x,x') = g-^/\x)Q{x,x')g~^/\x') 



g ^^'^{x)Q{x, X ) 



(6.10) 



Substituting (U),(|6;|) and (|O0|) into we find 



m'+(^--]R 
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53 - ^Rjy'^dpQ 



1 1 1 \ 

+^R\xpR^\'^ y"y^y^dsQ 
+ [^R\''0;,s + ^R\xpRVs) y'^yWd.d^Q = -5(y), (e.ii) 

where are the coordinates of the point x and — {d/dy°'). 

We have retained only terms with coefficients involving four derivatives of the 
metric. These contributions give the ultraviolet divergences that arise in the course 
of renormalization. 

In normal coordinates with origin at x' , ^{x, x') is a function of y and x 
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Q{x,x') = Q{y,x') (6.12) 



where y belongs to the small region around x' . In this way the equation ( |6.6| ) may be 
solved recursively. Namely, we will introduce the momentum space associated with 
the point x' {y = 0) by making the n-dimensional Fourier transformation: 



^ix,x)= I j^-y^^ik) exp{iky) (6.13) 



where ky = kay"^ = Tj^^kay/s and expanding '^{k) in a series: 



Q{k) = %ik) + Qiik) + Q2ik) + ... (6.14) 



or 



9, 



{x,x') = J j^^Qi{k)exp{iky), i = 0,1,2,... (6.15) 



where we will assume that the coefficients Q{k) have geometrical coefficients involving 
i derivatives of the mertic. 

On dimensional grounds, '^i{k) are the order k~^'^~^^\ so that ( |6.14D is an asymp- 
totic expansion of Q{k) in large k (small y). 



Inserting ( |6.13|) into (|6.11|) we get that the lowest order solution (~ Oiy"^)) is 



%{k) = {k^ + mY (6.16) 
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and 



Qi{k) = 



The function '^2{k) (~ 0{y^)) may be found from 



{r^d^dp-m')%-{^--)R% 



Using ( |6.13| ) one can get that the last two terms of ( |6.18| ) cancel each other. 
In another way we can consider that Q^o is Lorentz invariant of the form 



Then, inserting ( |6.19|) into (|6.18| ) we get 



- ^Rjy''d(,%iy) + ^Rj,%''y''dpds%iy) ^ 0, 



and 
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Therefore 



(6.22) 



The Lorentz invariance of ^oiu) leads to further simphfications of ( |6.11|) . Namely, 
the contributions 



-^Rj;^ + IRcT^) y^y^dpC^iy) + ^Rx\^,^yVy%d^%{y) = (6.23) 



are eliminated. 



In the same way 



— —R a-.B-y + —Ri 



20 10 



+ ( ^R^'^p^.s + ^R\xpR\'')j y'^y^y'^y'd^d^^.iy) ^ 



(6.24) 



and we have the following equation for '^{y) to the fourth order in derivatives of the 



metric: 



%y) 
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1 

30 



60 



1 

40 



-5{y) (6.25) 



Substitution of '^2{k) ( |6.22| ) instead of ^^{k) in the identity ( |6.20| ) does not change 
it, thus we can suggest that ^2{}i) is Lorentz invariant too and it can be wtitten as 
^2(?/) ~ iy'^ya)'^- The equation (|6.11|) is simphfied to 



6 



6 



+ 



'1/2) [i - R;aP - -^Ra^Rp-r + ^Ra'^ -y^ RpS 



d^'d'^Qik) = 1 (6.26) 



where 



5°53(A;) = dQ{k)/dk 
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Making a further recurrent process with this equation we get 



and 



where 



%{k) = (6.27) 



1 ^2 



+ a^p{k^ + mY^d^'d^ie + (6.28) 



aaP = (1/2) (e - I) R;aP + 



n 5. 

60' 



Now we can write the equation for the Green 's function in a convenient form. 
Let us introduce useful equations: 

{k^ + mY^d'^ik^ + m^)-! = (l/2)a"(A;2 + m^)-^ 
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and 



2 I _2\-3 



(6.30) 



Using ( |6.30| ) we may write Q{k) as 



+ {l/3)ao,pd''d^{k^ + m 



2\-2 



+ 



(A;2 + ^2)-3 



(6.31) 



Inserting the last equation into (|6.13| ) we get 



d^k 

'^{x\y) = I j^^i^ exp{iky) 



1 + 7i(x',?/) 



d 



dm? 



d 



dm? 



(6.32) 



where, to the fourth order in derivatives of the metric, the coefficients are: 
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li{x,y) = {^-^^R+ (1/2) - e ) Ray'' - (l/3)a„^y V- 



72(x',i/) = (1/2) Q - - (l/3)a^ (6.33) 



The expression for the Green's function ( |6.10D will be then | |Bunch fc Parker 1979 



/ r dT'k 

G{x ,y) = g'^^'^iy) J j^^exp{iky)x 



2 / a \:) 



E7,(x,y)(-^j {e + mY' (6.34) 



X 

j=0 



The equation ( |6.34|) is very important for further calculations. 



6.2 The Green's function and the 
Schwinger-DeWitt method 

We have found already that the geometrical quantities enter directly into the 
structure of Green's functions for arbitrary fields through the covariant d'Alembertian 
operator and non-minimal connection. Now we will treat the equation for the Green's 
function with the Schwinger-DeWitt method developed for flat space-time in the 
previous section. 

Let us multiply the equation ( |6.6|) on the left side by g^^^^\x) and on the right by 
^(1/4)^^'^^ and introducing 
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^{x, x) = g^^/'^\x)G{x, x')g^^/''\x') (6.35) 

we will rewrite this equation in the form 

(-df.d" + + ^r) ^{x, x') = S{x, x), (6.36) 

where S{x,x') = g~^^'^{x)S{x — x) is scalar with respect to general coordinate trans- 
formation, and product of 5 -functon is 

(1, 5{x, x')) = j (fx5{x, x) = 1 (6.37) 

For simplicity we put ^ = and get; 

(-U^ + m^))'^{x,x) ^ 5{x,x), (6.38) 

or, in operator form 

F-^s = 1 (6.39) 

where F is matrix operator. 

Let us introduce the representation for $5 in the form 

oo oo 

'^s{x,x ; s) — i < x\ j dsexp{isF)\x >— i j dsf{x,x ;s)exp{—im^s) (6.40) 
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We can get from ( |6.38| ), that the function f{x,x'; s) is the solution of the equation 



d 

n^f{x,x]s) = i—f{x,x']s) (6.41) 



where information about space-time structure is included in the d'Alembertian. 

We may turn this equation into an elliptic one by rewriting the equation for 
amplitude with the replacement x^^^ = ix^^'^ and s = it. We will have 



^^f{x,x'-t) = nj{x,x';t) (6.42) 



One can write the solution of the ( |6.4(]| ) as a simple expansion of the solution for fiat 
space-time. This solution is 

f{x,x';s) = (47rs)-'^/^exp(-|x-x'|V4s) (6.43) 

Returning to the initial variables we get 

f{x,x';s) = (47rzs)""/^exp(-|a; - x'|V4is) (6.44) 

In curved space-time we may expand this solution to a local asymptotic expansion 
(for X ~ x' and s ~ 0) 



f{x,x';s)r^{47cisY "/^^ exp(-cr(x, x')/2is) ^ 7j(x, x')(is)^' (6.45) 

j=0 
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where a{x,x') is the so-called geodesic interval (half square of geodesic distance be- 
tween points X and x'). In particular 



/(x,x;s)~(47r^s)(-"/2)^7,(x)(^s) 

j=0 



(6.46) 



The explicit form of fj{x,x') can be calculated recursively Prown 1979| 



DeWitt 19751 , Fulling 1989 1 



In the limit of coincidence x ^ x one finds: 



7o(x') = l; 7i(x') = (l/6-0i?; 



and 



j,{x) = {l/6-0'R'-{mai 



(6.47) 



From the equations (|6.40|) , (|6.45| ) and (|6.35| ) we get the explicit expression for Green's 
function in the Schwinger-DeWitt representation: 



GsDix,x 



iA^/'^{x, x 
(47r)"/2 



exp (^—ism^ — a{x, x )/2isj ^ x ){isy 

j=0 



(6.48) 



where 
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A(x, X ) = —g{x) ^^'^det[dadf3a{x,x )]g{x) ^ 



-1/2 



(6.49) 



is the Van Vleck determinant fin normal coordinates about x this determinant is 



reduced to g-^/\y) (|04D). 



6.3 Connection between the two methods 



Let us put 



(/c^ + m^) ^ = j ids ex'p[—is{k'^ + m?^ 



(6.50) 



Then integration over the momentum in ( |6.34| ) leads to 



I exp[— is(fc^ + w?) + iky] = i{ATris)"'''^ exp(— zsm^ — a/2is) (6.51) 



(27r) 



The resulting equation for the Green's function will be 



X exp 



—is{k + m ) + iky F{x ,y,is 



(6.52) 



where 



F{x',y,is) = J2lj{x\y){ 

j=0 
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Comparision with the expression of the Green's function in the form of ( |6.48D gives 
A(x, x') = g^^^'^iy), and 

,y) = -fj{x,x') (6.53) 



In this chapter we considered two methods for computing the Green's function of a 
scalar field in curved space-time and we prepared the basis for future finite tempera- 
ture calculations. 
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Chapter 7 



FINITE TEMPERATURE 
BOSONS 

IN CURVED SPACE-TIME 

In the previous section we developed a mathematical formalism which is a convenient 
tool for the description of thermal Bose gas in curved space-time. In this section we 
can consider the ensemble of bosons interacting with gravity at finite temperature. 
Let a total system ("matter and gravity") be described by the action 

The gravitational action is 



Sg = j d'^x^Jg{x)Lg 
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(7.2) 



with Lagrangian 



and the action for a matter field is 



To write the action for the matter field we use a generating functional Z[J] 



Z[0] = J Dipcxp i^-{t/2) J d^x^/^Mx){-n, + m^ + ^RMx: 



oc DetG{x, x'Y^^ exp {-{i/2){J{x), G{x, y)J{y))) 



Then the functional W[Q\ = — ilnZ[0] will be 



W[0] = -{i/2)lnDetG{x,x' 



— (i/2) j ids{is)~^trf{x,x ,is) exp^—im^s) 



We may connect the Green's function and heat kernel with the equation 



Gsd{x,x ) = j idsf{x,x , is) exp(— im^s) 
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and write ( [7.6[ ) in the following form 



W[0] = -{i/2) / d^x^gix) / dmhiGsD{x,x') {71 



From the equation ( [7.8|) we find an important expression for the effective Lagrangian 
of matter field 



Leff,m. = (— ^2) / d'^xyg{x) / dm'^tTG{x,x ) (7.9) 



Then the effective action will be: 



Seff = / d^x^/g{x)Leff{x) 



d\/g{x) Lg{x) + (-z/2) / dmhiG{x,x') (7.10) 



In order to apply the usual formalism of finite temperature quantum field theory, 
we will assume that the space-time manifold M4) is a static manifold with topology 
X M3) where refers to time coordinate and M3) is the spatial, three dimensional 



section of M4). We will choose M3) without boundaries [[Kennedy et al. 1980|| , then 
no surface terms will appar in the induced action. 

The heat kernel may be expressed as the sum of zero-temperature images 
Denardo fc Spalucci 1983], Powker fc Kennedy 1978], Powker fc Critchley 1977[ 



119 



/(a;,x,zs)= 2^ . /3) X,x,«s 7.11 



where the sum goes from periodic restrictions and f3){x,x,is) is the solution of the 
three-dimensional equation 



''¥s ~ °" ^ H /3)(5,5,^5) = (7.12) 



This solution is written in the form of a series 



fs)ix,x,ts) = J2^,iR)itsy (7.13) 

n=0 



Then from the equation (|7.7|) in the limit of coincidence (x — > x ) 

oo 

limG(x, x')= hm idsf{x,x\is)exp{—im'^s}. ('''•14) 

x—^x' X— >a:' J 


we get the finite temperature Green's function in the Schwinger-DeWitt representa- 
tion: 

hm G^5£)(x, X ) 



oo oo 

6 



/"OO 

Y^^ji^') J ids{isy~^/^exp{-ism'^ -n'^P^Ais) (7.15) 



(4vr)3/^±^oo,eo '^^ ^Vo 
Selecting the temperature independent part {n = 0) we find 
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lim G^sd{x,x') =Gsd{x',x) + G,>{(3) (7.16) 



where the finite temperature contribution G^> (/3) is 



^ oo oo 



G.'(/5) = 7T-T^2EE7,(X')X 
J j=On=l 



oo 



X / ids{isy'^exp{-ism-n^(3y4is) (7.17) 







and Gsd{x ix) is the hmit [x x ) of the Green's function in the Schwinger-DeWitt 
representation (|6.48|) . 



Summation in ( [7.17| ) may be done with ( |3.27| ) and ( p.31 
In the result we find 



GAP) = 77^ Y.ElA^')iM^^y-'K,.si(3mn) (7.18) 

l^^J n=l 9=1 



Total Lagrangian may be written in the form; 



oo \ oo 



L,jf{f3) = \Lg{x) + (-V2) J dmhrGsDix,x) + (-z/2) J dm^GAP) (7.19) 

The first two terms give the same geometric structure, and after renormalizations we 
will have the Lagrangian L^, the third one is temperature contribution f{(3). 

L,ffi/3) = L,-fi(3) (7.20) 
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where 



Lg = Lg{x) + (-■ i/2) / dm^iiGsD{x,x 



(7.21) 



and 



f{(3) = {i/2) / dm'GAf^) 



(7.22) 



The last expression can be written as series 



j=0 



(7.23) 



where 



2 CXD 



^o(/5H = -TT-^ Y.(l/n^)K,{f5mn) 

P n=l 



2m °° 



2 °° 

62 (/?m) = -— ^ Ko{Pmn) 



n=l 



(7.24) 



Using the results of Appendix: ( pCX11.28D , ( pCXlLSOj ), ( pCX11.32D , we may rewrite 
(|7. 23|) in the form of series [Kulikov & Pronin 1993] 



(7.25) 



122 



where jj{R) are ( |6.33[ ), coefficient bQ{/3m) is 



fco = (l//5) / ^ln(l-exp(-/56)); 



(2vr) 



and 



hj{(3m) 



d 



dm? 



bo{/3m) 



(7.26) 



The expression bQ{/3m) in ( |7.25| ) is the density of Helmhohz free energy in flat space 
and the following ones are created by corrections which are connected with the inter- 
action of the heat bosons with gravity, so the equation ( [7.23| ) describes the Helmholtz 
free energy of a Bose gas in curved space time. 
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Chapter 8 



FERMI FIELDS 



IN CURVED SPACE-TIME 



8.1 Momentum-space representation 

for the Green's function of a fermion 

Now we will get the Green's function of fermions in curved space-time to apply 
the above results for computation of the effective action of the system "matter field 
+ gravitational field" . 

As we know already, the fermionic action with Lagrangian ( 5.40| ) is written as 



= (z/2) / d^x^/^)^ix) (t^V^,, + m) iIj{x) (8.1) 



The generating functional then is 
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Z[v, v] = Di/jDi/j exp iS^ + i{r],^lj) + i (^ip, r] 



oc {DetGp) ^^"^ exp \—i {ri{x), Sf{x, y)"!])] 



(8.2) 



where the scalar product includes the square root of the metric. 

Connection between the bi-spinor Gp and Fermionic Green's function Sp is 
Birrell & Davies 198211 : 



Sf{x, y) = in'^D^ + m)GF{x, y) 



(8.3) 



and the Green's function satisfies the equation 



[i^D^ + m)SF{x, y) = -g{x)-^'^5{x - y)l 



To solve this equation with momentum space methods we introduce Riemann normal 



coordinates [Petrov 1969 1 and write the spin connection in the form 



Panangaden 1981 



T,{y) = (l/16)[7„,7/3]^"V2/'^ + 0(y^) 



and the vierbein field as 



K{y) = {i/Q)t''R,upcy^y^ + 0{y'^ 



(8.6) 
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The expression for 7^(x') may be written as 



The spinor derivative appearing in Dirac's equation is written as 



(l/16)7l7a,7/3]i?"V2/' (8-8) 



The Fourier transform of S{x\ y) is 



— — exp(z%)5(fc) (8.9) 



(27r) 



Let 



5(fc) = 5o(fc) + Si{k) + 52(fc) + ... (8.10) 

be an asymptotic representation of S{k) for large k. 

The values Si{k) are asymptotic variables of the order k~^^~^^K They may be found 
with recursion procedure from the equation 



126 



-(l/16)7'^[7.,7/3]^"Vl/'^ + -1 S{x,y) = 5{y) 



(8.11) 



In this case the momentum space representation of a propagator of a fermion will be 



dP'k 

S{x\y)= I j^exp{iky) 



^2 _|_ jjj2 (J^2 _|_ 1^^2-12 



+ ... 



(8.12) 



The momentum space solution of the equation for the bi-spinor 

(□^ + l/4i? - m^) Gi.(x, x') = -g{x)-^/^5{x - x')i (8.13) 

where = D^^D^^x is the covariant d'Alembertian of spinor field, may be obtained 
with the same momentum space methods as in chapter VI punch fc Parker 1979| . 



The result can be written in the form of the following expression 



Gf{x, X ) = Gf{x , y) 



r d^k ^ , / d 



{k' + m')-^ (8.14) 



where the geometrical coefficients aj{x\ y) in the limit of coincidence {x ^ x) are 
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aoix\y) = i; 



a,{x,y) = {l/12)R-l; 



and 



a2{x\y) = (-(l/120)i?^''^ + (l/288)i?2 



(l/180)i?^,i?^^ + il/180)R^,^rR'"'n ■ 1 



+ (l/48)S[„,^]S[,,5]i?"^^«i?7A' 



(8.15) 



where = (l/4)[7a, 7/3] 

8.2 The bi-spinor function in the 

Schwinger-DeWitt representation 

In analogy with the scalar field we can rewrite equation ( ^.141 ) in the Schwinger- 
DeWitt representation. Since 




(8.16) 
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the equation ( ^.14| ) will be 



Gf{x, X ) = — 77-^7j^ — / exp 



(4W2) 7 (2s)"/2 



+ m2) F{x,x ■,s) (8.17) 



where 



(8.18) 



n=0 



and coefficients d„(x, x') are defined by ( ^.15| ) 

As in the scalar case the determinant A'-^^'^\x, x') is defined by the equation 



A^'/'\x,x)=g-'/'iy) 



(8.19) 
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Chapter 9 



FINITE TEMPERATURE 
FERMIONS 

IN CURVED SPACE-TIME 

9.1 The Helmholtz free energy of a Fermi gas 
in curved space-time 

After we constructed the free energy for a thermal scalar field we can consider a 
thermal fermi field. Let the total Lagrangian of the system of "gravity + fermionic 
matter" be 

Stot = Sg + Sra (9-1) 

where Sg is ( |7.2| ) and 
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^m = J d'^x\l g{x)Leff^^ (9.2) 



To write the action for a spinor field we will use the functional method for calculation 
of the generating functional from ( |9.2|) . Making the same procedure as in ( [7.6| )-( [r8D , 
write 

W[Q\ =ilnZ[0] = {i/2) \nDetGF{x,x') 



ids{is) ^tTf{x, X , is) exp(— -im^s) 



{i/2) j d^x r dmhrGF{x,x') (9.3) 



where the Green's function is expressed by 



Gf{x,x)= / ids{is)~^ f{x,x ,is) exp{—im'^s) (9.4) 
Jo 

and f{x,x',is) is the heat kernel. 

The kernel f{x,x',is) is (in the limit (x — > x')) the sum of zero-images 
Uenardo fc Spalucci 1983"| antiperiodic in the imaginary time 



ff ' M exp[-/?^(n-l/2)V4z.] ^^ 

f{x,x,is)= 2^ . f3){x,x,is) (9.5) 

n=-oo (47r«s)V^ 
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where an equation for fs){x, x\is) is 



'^^-a3)-(l/4)i?]/3)(i,5,^s) = 



(9.6) 



and is the covariant d'Alembertian on M3). 
The solution of (WM) is the series 



f3){x,x,is) = ^dj(a;')(zs)^ 

j=Q 



(9.7) 



where the coefficients aj{x') are determined by ( ^.151) . 

From the equation ( |9.4| ) we get the finite temperature Green's function 



hm^ Gpsij{x, x) = Gf,sd{x ,x) + G'f(/5) 



(9.^ 



where Gf,sd{'^\x ) is ( |8.17| ) in the hmit {x = x ) and 



X j ids{isy-^^'^ exp [-iW - {(3^/Ais){n + l/2f 




(9.9) 



Summation with respect to (|3.27|) with z = 1/2 gives 



(9.10) 
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(9.ii: 



After integration over the proper time (s) ( |3.31|) we get the following expression for 
the finite temperature contribution in the Green's function of a fermion 



= 7X7^ E E a,{R){-iy{Pn/^my-'K,.s{(3mn) (9.12) 

l^^J j=On=l 



The total action of the system will be 



Uff{P) = L,{R)-MP,R) 



(9.13) 



where L{R)g is the temperature independent Lagrangian 



L{R)g = Lg + (i/2)tr / dm?GF,sD{.x,x ) 



(9.14) 



and the finite temperature contribution is expressed in the form of a series: 



fF{f3,R) = (-V2)tr / dm^GF{l3) = Y.a,{R)f,{Pm) 



(9.15) 



with coefficients 



aj{R) = (l/2s)traj(i?), 



(9.16) 
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and 



2s °° 

71=1 



(9.17) 



The finite, temperature dependent contribution R) represents the density of 
Helmholtz free energy in curved space-time. Using the integral representation for 
the series of modified Bessel functions ( pCXll. lO] ) , ( pCXll.T^ ) , ( pCXll. 141 ) one can write 
it as [Kulikov & Pronin 1995]: 



MP, R) = foiPm) + a,{R)h{Pm) + a2{R)h{M + • • (9-18) 

where the first term is the standard form of the Helmholtz free energy in Euclidean 
space 

/o(/3m) = -|/^0jln(l+e-*) (9.19) 

with energy of particle e = \JkP' + m?. 

The factor 2s = 4 reflects the existence of the four degrees of freedom present in 
the fermion field: particles and antiparticles, spin up and spin down. 
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The following terms are geometrical corrections of the Riemann space time struc- 
ture with respect to the Euclidean one with temperature coefficients in the form 

The method developed above does not allow us to compute the density of the grand 
thermodynamical potential. Therefore, in the following calculations we will use the 
local momentum space formalism as the most convenient for the construction of local 
thermodynamics. 
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Chapter 10 

THERMODYNAMICS OF 
GAUGE FIELDS 

In this chapter we will apply the formalism of gauge fields in curved space-time 
developed in chapter V to study the properties of thermal photon gas in an external 
gravitational field. 

10.1 The Green's function of photons 

As we know already from chapter IV the total Lagrangian for a vector field in 
Minkowski space-time is the sum of three contributions: 

Ltot — Lm + Lf + Lgh, (10-1) 
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where 



= -(l/4)F^,F^^ 



and 



Lf = -{l/2a){d,A^f 



It can be extended to curved space-time with the transformation 



where = + and the connection is defined by equation 



= (l/2)E„^/i"^(x) 



d 



h^nU{x) 



where the matrices T,ap are 

The strength tensor of electromagnetic field is 
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The variation of the action 



S = j d^x^{L + Lf + Lgh) (10.8) 



gives the equation for the vector field: 



V.F'^" + (l/a)V^(V,A'^) = (10.9) 



In the Feynman gauge (a = 1) this equation has the form: 



V^VM'^ - V^V^A^ + V^V^A^ = (10.10) 



and with the definition of the Riemann tensor ( |5.17| ) it is 



V^,V''A^' - R^'^A'' = {) (10.11) 

Based on this equation we write the equation for the Green's function 

V^VD^ix, x') - R>'^D\{x, x') = -g-^/\x)6{x - x')5^ (10.12) 

We will be interested in calculations of the Green's function in the limit (x ^ x) to 
find the effective action. 



We will rewrite the equation ( 10.12 ) in Riemann normal coordinates with origin 



at the point x' . For covenience one may define the Green's function Dl^(x,x') as 
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= g-^/\x)D'^{x,x') (10.13) 

where we used the fact g{x') = 1 

For our calculations we will use the Christoffel symbols which in the Riemann 
normal coordinates are 

= -(l/3)(i?^^^ + R%a^)y^ (10.14) 

where y'^ represents the coordinates of the point x and the Riemann tensor is evaluated 
at x . The expansion of equation ( |10.12|) to the second derivative of the metric gives 



r^d^dpD^^iy) + {l/6)RD'^{y) 

- {4/3)Rmy) - {i/mly'dxD^M 

+ {l/m%\y'y'd.dpD^:{y) - {2/3)R^,\y'd^DJ{y) 

+ {2/3)R^\^yW^D^^{y) = -5{y)5'^ (10.15) 

where da = d/dy^. 

The momentum space approximation is defined by introducing the quantity Dl^{k) 
defined as 
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D'^{x,x') = D'^{x\y) = j ^D^f(fc)exp[z%] (10.16) 

This quantity is assumed to have the expansion 

D^,{k) = Di;^,{k) + D^,^,{k) + Dl^^^ik) + ... (10.17) 

where D^^{k) have a geometric coefficients involving i derivatives of the metric. On 
dimensional grounds D^^(k) must be of order fc"*^^"*"*) so the equation ( 10.15|) is an 
asymptotic expansion in large k. 

As we can see from (|10.15|) the first term in ( p,0.17| ) is 



D^^M=S'^/k', (10.18) 



The second one is 



DW) = 0, (10.19) 



The third one is computed from 
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+ {2/3)R^\^yWaDUy) + ... = 



(10.20) 



To simplify calculations note that y° — > {—id/dk°'), and integrate by parts to find 
that 



(27r 



(27r) 



X 



{{l/Q)R5^r - {2/?>)R'X} jk^ - {A/3)R^'',^k^kp5';/k^] (10.21) 



The final expression for the photon propagator will be then: 



7|^exp[iA;y]x 



(27r)^ 



X 



+ {(l/6)i?5^^ - (2/3)i?(:(5^} /k^ - {A/?,)R^^,^k''kp5';/k^] (10.22) 



The Green's function of the ghost fields obeys the same equation as a scalar field, 
therefore one can use this equation for our further calculations. 



10.2 The thermodynamic potential of a photon gas 



As in the case of fiat space-time, one must carry out the calculations of the free energy 
of a photon gas together with the ghost contributions. 
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The generating functional of an abelian vector field in curved space-time is written 



as 



X exp 



Z[J^,r7,r/] = J DA^DrjDrj 



(10.23) 



Integration over the fields with zero sources leads to the following result for the 
logarithm of the generating functional 

In Z[0] = In Det{G) - (1/2) In Det{G^iv) (10.24) 



where the Green's functions are defined by ( |10.22| ) and ( |6.34| )p| 



From this equation and from the definition of the free energy directly follows the 
expression for density of free energy 

oo 

fph{P,R) = {-i/ 2) liin J dm^[2trG{p,x-x')-tiGf,^{p,x-x')} (10.25) 
The final result may be found after putting a mass parameter into the expressions 



for the propagators of the photon and ghost fields and, after calculation of ( |10.25| ) 
setting the mass equal to zero. 



^To find the Feynman propagator of tlie gliost field from tlie boson propagator we liave to put 
the mass of boson m — after calculation. 
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After making this procedure we will have 



U,{P,R) = Y.g,{R) f-^ytrln(^^ + e2) (10.26) 



J 



where the symbol tr means 



tr= 2^ 

n=— oo 



(2vr) 



3 



and gj{R) are geometric coefficients. 

After summation in (|10.26|) we will have the final result in the form 



~{l/Q){R-2R^j;)[e{exp[l3e] - 1)]-'} (10.27) 



This expression is the density of the Helmholtz free energy of a photon gas in an 
external gravitational field. 



10.3 Internal energy and heat capacity of a photon 



gas 



Thermodynamic properties of the photon gas have been studied very well in fiat 
space-time [[Landau fc Lifshitz 1959|, Puang 1963||. The formalism developed in this 



section allows us to get the properties of photon gas in curved space-time. 
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The density of energy of a photon gas is 



2e 



exp[/9e] — 1 



iimR-2WX^)-j—-^-^\ (10.28) 



The first term in ( |10.28| ) corresponds to the results of statistical mechanics, and the 



second one is the curved space-time correction to the energy of photon gas. Integrating 
over the momentum 



(fk 1 1 



(27r)3e(exp[/?e] - 1) 12(3 
we get the following approximation to 0{R^) 



(10.29) 



u{R) = aT^ + (l/72)(i? - 2R''^5l)T^ + ... (10.30) 



where a = (vr^/lS). 

The heat capacity will be 



c{R) = ^u{R) = AaT' + (l/36)(i? - 2R^J;)T' + ... (10.31) 



As was shown in [[Altaie fc Dowker 1978|| the effects of curvature for Einstein space 



lead to the results for the Planck black body expression which are equivalent to the 
results (p^:30|) . 
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Chapter 11 



RENORMALIZATIONS IN 



LOCAL STATISTICAL MECHANICS 



11.1 Divergences of finite temperature 
field models 

As one can show by calculations, the contributions Gsd{x,x) in the expressions 
( |7. 161 ) and ( |9.8D are divergent [[Birrell fc Davies 1982|| . To have a clear picture of the 



model under consideration we need to eliminate these divergent contributions. The 
direct way to eliminate the divergences is to combine the gravitational Lagrangian 
and the divergent parts of the matter Lagrangian. 
1) Renormalization of a hose field. 



As was shown in ( [7.21|) the effective Lagrangian for gravitational field can be 
written as 
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Lg = Lg{x) + {-i/2) / iiGsD{x,x 



11.1) 



Inserting (|6.48|) into the last term of the expression we find 



-z/2) j iiGsD{x,x) 



2(47r) 



1 ~ 7 



11.2) 



We will use the procedure of dimensional regularization |ft Hooft fc Veltman 1972] to 



select divergent terms of the expression ( |11.2| ). It gives 



(-z/2) / iiGsD{x,x) 



1 oo » 

-(47r)-"/2^7,(a:)/^rf.( 



zs)(^'~^-"/2)exp(-zsm2) 



1 / 777 \ OO 

|(4vr)-"/^ (m) E 7,(x)m^-2^r (j - n/2) 



;il.3) 



where n is a dimension of the space-time and M is an arbitrary mass scale. This 
parameter is introduced to preserve the dimension of the Lagrangian to [L]''^ for the 
dimensions n 7^ 4. 

The functions T[z) have the poles at n = 4: 
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n-'i 



2 J nin - 2) \A-n 



7) + 0{n - 4) 



r 1 



11.4) 



Therefore the first three terms of the expression ( p.l.3| ) are divergent. 
Selecting divergent parts and using the logarithmic expression for 



n— 4 



1 , m 



l + -ln^ + 0((n-4)2) 



;il.5) 



write the divergent contributions in the effective Lagrangian L„ in the form 



{-1/2) / trGsD(x,x') 



-(4vr) 



-n/2 



1 1 

+ 



n-4 2 



j + In 



m 



X 



4m 



— 4) n — 2 



m^7i(x') + 72(x') + 



11.6) 



Combining (|11.6|) with the gravitational Lagrangian ( [7.3|) we can redetermine the 
constants of the gravitational Lagrangian (|11.1| ) as 



1a 1a 



1 1 m 



ttGr ^ SttGo (47r)2(n-4) 2 



11.7) 
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1 1 1 1 2^1 . 

+ 7T^ in-^ 



IGttGb^ levrGo (47r)2n-4 V6 



180 (47r)2n-4 



7fl = 7o - ./s^ ^ . (11-8) 
180 (47r)"^ — 4 



where Go, Aq, oq, /?o, 7o are bare constants and Gr, A/j, a^j, are the physical 

(finite) constants. 

Since all divergences of the matter field can be included in Lg, only the finite 
temperature contributions observed in ( [r.20D remain. 



2) Renormalization of a fermi field. 

As follows from ( |9.14|) the effective Lagrangian of the gravitational field is 



L, = L,-^(47r)-"/2^trd,(i?) 



2 



X 



X / irfs(is)^~"/^"^exp(-iW) (11.9) 



After dimensional regularization procedure (as in scalar field case), we get 



1 . 1 1 2m^ 

-A/? = TT^Ao + 



SttGr " 871^0 (47r)2n-4 
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1 _ 1 1 1 

IQ-kGr ~ levrGo " (47r)2 (n - 4) T 



1 1 1 



(47r)2 (n -4) 144 



(47r)2 {n - 4) 90 



Thus, the Lagrangian ( p.l3| ) does not include the divergencies in the finite tempera- 
ture loop approximation over the fermi fields. 

The method of dimensional regularization developed here is not the only one for 
applications to field models in background curved space-time. The (^-function method 
[ [Hawking 1977|| , [purt 1983|| and the method of covariant geodesic point separation 



also lead to a solution of the problem of regularization of effective Lagrangians of mat- 
ter fields and calculation of Energy-momentum tensor anomalies ||Christensen 197^ . 
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Chapter 12 



LOCAL QUANTUM STATISTICS 



AND 



THERMODYNAMICS OF BOSE GAS 

In this chapter we will continue to develop the methods of local statistical thermo- 
dynamics [[Kulikov fc Pronin 1995|| with application to the ideal bose gas in curved 



space-time. In chapter VII working with the Schwinger-DeWitt proper time formal- 
ism we found that the density of Helmholtz free energy in curved space-time may be 
found in the form of a series ( [7 .251 ) including the geometric structure of space-time. 
However, as in the previous section, working with this formalism we could have cer- 
tain difficulties with the introduction of the chemical potential and the computing of 
the density of the grand thermodynamical potential. 

To avoid such problems which will appear on the way of construction of local ther- 
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modynamics with the Schwinger DeWitt formahsm, we turn to the local momentum 
space method [ Bunch fc Parker 1979 |, [ Panangaden 198T | in quantum field theory in 
an arbitrary curved space time and the imaginary time formalism to introduce the 
temperature jPolan fc Jackiw 1974|| , [[Weinberg 1974 |. 

For the construction of the local quantum statistics and thermodynamics of a 
bose gas, we will use the connection between the partition function of ideal quantum 
systems and finite temperature Green's functions, which may be found by the local 
momentum space method. 



12.1 Density of grand 

t her mo dynamical potential 

In this section we will analyze the scalar field model with a conserved charge. The 
Lagrangian of the model is 

Sm = -(1/2) J d^x^)^*{x){-U, + m^ + ^R)^{x) (12.1) 

where $ = (0i, 02) is a dublet of the real fields. 
The action written in terms of real fields will be 

Sm = -(1/2) J d^x^^)(P\x){-n^ + m^ + iR)r{x) (12.2) 

The total action of the system "matter + gravity" is 
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Stot — 'S'g + Sm 



(12.3) 



Now we can write the effective action at finite temperature in analogy with ( [7.20| ) as 



;i2.4) 



where Lg is (|7.21|) and uj{P,fi,R) is the density of grand thermodynamic potential. 

The result ( |12.4|) may be obtained with the momentum space representation for 
the Green's function of a boson ( |6.34| ) fKulikov fc Pronin 1987| . 

In the momentum space representation, the expression for Lg// is split into two 
parts 



Le// = — (i/2) J dm'^tTG{x, X ) — u{/3, n, R) 



;i2.5) 



The potential a;(/3, R) is 



2 oo 



dm2^^7,(a;,x') 



d 



j=0 n=0 



dm? 



X 



X / (0K^ + a-^ 



(l/2)E7.(a:') f- 

.1=0 \ 



d 



tr ln(ct;„^ + e^^ 



dm? 



;i2.6) 



and coincides with the Helmholtz free energy ( 7.22|) and ujn = 2Tm/ (3 
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The symbol tr in ( |12.6| ) is determined as 



tr(...) = j:[ 



For introducing the chemical potential, we will change the Matsubara frequencies 
Un ujn + fJ' and thermodynamic potential will be /x, R). 

Since both positive and negative frequencies are summed, we will have 



trln(u;„^ + e^) trln[(cj„ + /i)^ + e^^ 



tr |ln ujn^ + (e — I^Y 



+ In 



;i2.7) 



After summation in ( |12.7|) we get 



(12.8) 



where 



UJ. 



.{f3, /X, R) = (1/(3) E 7,(x') j^-^y ln(l - exp[-/5(e - /i)]) (12.9) 



and 



u+{(3, /i, R) = (1//5) ^ ^,{x') (^-^y ln(l - ^ exp[-/5(e + /.)]) (12.10) 
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So, with accordance to ( p.2.8| ) the density of grand thermodynamic potential is the 



series 



= j2l',{x")b,{Pm,z) (12.11) 



where 



ho{(3m,z) = {1/(3) ln(l-zexp(-/3e)); 



d 



1 



b,{(3m,z) = [-g-^j bo{(3m,z) (12.12) 



and the fugacity is ^ = exp(/5/i). The geometrical coefficients 7j(-R) of the equation 
( |12.12|) have the form (|6.47|) . Renormalizations in the total Lagrangian (|12.4| ) are the 



same as in chapter XI. 

12.2 Statistics and thermodynamics of bose gas 

We find the bose distribution function as the derivative of the grand thermodynamical 
potential 

For occupation numbers with momentum k we get 
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n 



k [z-^el^'i - 1 



;i2.i3) 



where the function -B(/3, R) is described by the formula 



B{P,R)= l + 7i(i?) 



A. 

2er 



1 - ri - ze-^'^^)-^ 



+ 



;i2.i4) 



The function -B(/3, i?) depends on the curvature, temperature and energy of the boson. 

Studying the thermodynamical properties of Bose gases we will start with the 
equation 



UJ 



(12.15) 



In the non-relativistic limit for particle energy e = k"^ /2m we get from (|12.15|) the 
equation 



uj{(3,n,R) = J2lj{x')95/2{ 

j=0 



d 



dm? 



X 



(12.16) 



where A = (27r/mT)-'^/^ is a wavelength of the particle, and the function (75/2 (-z) has 
the following form 



00 



1=1 



15/2 



-(4/v^) J dxx^ In (1 — zexp( 



—X 



;i2.i7) 
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The average number of particles in a certain momentum state k is obtained as the 
derivative 



d 

< Uk >= -—uj{p,fj.,R) 
a/i 



j=0 



dm? 



z ^ exp(/5e) — 1 



(12.1^ 



The density of particles is 



n = 



l-7i(i?)(3/4m2)-72(i?)(3/W) g^,2{z) + (12.19) 



where the new function (73/2 (-z) is 



9Z/2K.Z) 



and no = zjiX — z) is the average number of particles with zero momentum. The 
functions (73/2(2;) and g5/2{z) are special cases of a more general class of functions 



00 
k=l 



12.20) 



In a more simple form the equation( p.2.19| ) may be written as 



{n - no)A^ = g^i2{z,R) 



(12.21) 



156 



where 



g^/2{z,R)= l-a^ + ... g^/2{z) (12.22) 



is a function which depends on curvature, and a is a numerical parameter. 

The equation (|12.21D connects four values: fugacity, temperature, density of the 



particles and curvature. 

We can solve it graphically and get the dependence of the (effective) chemical 
potential on curvature, temperature and density. The function (73/2 (-2, -R) varies with 



the curvature R as shown in Fig. I-l. The graphical solution of the equation (|12.21 
is presented in Fig. 1-2. In this Figure X-projection of the point of intersection of 
the curve ^'3/2 (-2, R) + A'^ng and of the line A^n shows the dependence of the fugacity 
on curvature. As follows from the graphical picture, for positive curvature i? < the 
fugacity z{R) > Zq (or fi{R) > yUo), for i? > we have z{R) < zq (or yu(-R) > yUo)- The 
fugacity zq (/iQ-chemical potential) corresponds to statistics in Euclidean space. The 
behaviour of the effective chemical potential is shown in Fig 1-3. 



12.3 Bose-Einstein condensation 

At low temperature there is a significal number of particles in the ground state, which 
can be expressed by the equation 

n^ = n-\-^g^,2{z,R) (12.23) 
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With the rising of the temperature the average number of particles no with zero 
momentum is lowered and for temperatures T > Tc it becomes zero. The temperature 
Tc is the critical temperature of a Bose condensation. The critical temperature may 
be found from the equation with {z = 1 and uq = 0) 



nXj^ = g3/2{l,R) 



12.24) 



The solution of this equation is [[Kulikov fc Pronin 1995 



2m? 



+ ... 



;i2.25) 



The temperature 



To = T,{R = 0) 



27r 



m 



n 



C(3/2) 



2/3 



12.26) 



where C(3/2) = 2.612... is the Riemann zeta function, n is the density of bosons and 
m is the mass of boson is degeneracy temperature (condensation temperature) in 
"flat" space without gravity. 
The ratio 



Tc{R) , R 



12m2 



[12.27) 



is not one, but depends on the curvature R of space-time. 
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As we can see, Tc{R) > Tq, for > and Tc{R) < Tq for R <0. The correction 

is small, therefore effects of curvature will be essential for quantum systems in strong 
gravitational fields. 
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Chapter 13 

LOCAL STATISTICS AND 

THERMODYNAMICS OF FERMI GAS 



The results of the previous chapter XII show the way to construct thermodynamic 
potentials of quantum systems with a variable number of particles. In this chapter 
we will develop this formalism for fermi systems. 

13.1 Grand thermodynamical potential and 

low temperature properties of Fermi gases 

An equivalent to the Schwinger-DeWitt representation, the momentum space rep- 
resentation of the bi-spinor Gf{x,x') has the following form ( ^.14| ): 




(13.1) 
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where 



Go{y) 



(27r)4 P + ^2 



(13.2) 



To introduce temperature, we will extend the time coordinate of the tangential 
space {y^} to the imaginary interval [0, —i(3\ and will consider fermionic field to be 
antiperiodic on that interval. Then, in imaginary time formalism, we can write 



f3 % 



{27r)*k^ + m'^ (3 J (27r 



d'^k g 



n=-oo ^n + ^ 



;i3.3) 



where ujn = nT{2n + 1), n = 0, ±1, ±2, . . . are Matsubara frequencies. To take 
into consideration the chemical potential we will shift the frequencies by cu^ ^ uJn + fi 
Morley 1978| , [Kapusta 1971 



Making the summation in ( |13.3 ) we will find Goly) in the limit of coincidence 
X = x' in the form 



1 



13 J (27r) n=-oo - (t^n + /i) + 



d^k \ I ^ ( e-fi 



e + n 



d^k f i (\ 



i I 



{2'nY\2e\2 exp /3 (e - /i) + 1 / 2^ \2 exp /5 (e + /x) + 1 



lim 



G+ {y) + G-p (y) 



(13.4) 
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This equation describes the temperature contributions for particles {/i) and antipar- 
ticles (— /i) separately. 

As a result one can find the expression for the finite temperature contribution to 
the Green's function for a fermion with non-zero chemical potential 

r d^k 2 / f) \^ , X -1 

'^''f^''') = / (^g«^(«) (S^j (1 + -"'T (13.5) 

Then the density of grand thermodynamical potential for fermions may be written as 

OO n 

uj{p,IJ,,R) = --tr J dm^Gp^x.ii) = ^aj{R)f){[3m;z), (13.6) 

where 



/o(/?m; ^) = / ^ ln(l + ze'^^) (13.7) 



and 



/, (/3m;z) = L^-y f,{(3m-z), (13.8) 



where 2; = e^^ is the fugacity and factor s = 2 (spin up, down). 

The coincidence of the finite temperature Schwinger-DeWitt and momentum space 
methods for calculation of the densities of themodynamical potentials is obvious for 
= 0. 
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Using the equations for thermodynamic potentials one can obtain some interesting 
properties of an ideal fermi gas in an external gravitational field : 

1) The Fermi distribution function of the gas in the gravitational field may be 
found from the expression 



riT 



duJk{f3,n,R) 



for occupation numbers with momentum k. 



It has the form 



nr. 



-le^^s + 1) 



(13.9) 



where the function F{P,R) is described by the formula 



F{f3, R)^ 1 + ai{R)^ [l - (1 + ze-^'^^)-' 



+ 



(13.10) 



and depends on curvature, temperature and energy of the fermion. 

2) We can estimate the dependence of the chemical potential on the curvature of 
space time in non-relativistic approximation. 

Let 



2m 



:i3.11) 
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then from 



n 



duj^P, /i, R) 



(13.12) 



we find the equation 



nA3(T) 



(13.13) 



where A = (27r/mT)^^^ is the thermal wave length of the particle, and 



"0 - 



3ai(i?) 3 a2{R^ 



4 m? 16 



;i3.14) 



is some function with respect to and n is the density of Femi gas. 
The function f-±{z) is 



/3/2(^) = E 



oo 



\n+l 



n=l 



n 



3/2 



dx- 



X 



^ -2 ^ exp(x^) + 1 



;i3.i5) 



The equation ( [13.13| ) may be solved with graphical methods. As we can see in Fig. 
1-4 the fugacity (chemical potential) depends on the curvature R of the space time. 

3) The explicit expression for the chemical potential at low temperatures and 
high densities {n\^ » 1), where quantum effects are essential, may be found by the 



calculations of ( 13.13 ) with the following representation of the function f^/2{.z) for 



large z [Huang 1963 
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f3/2{z) 



1/2 



+ 0(z 



-1\ 



(13.16) 



Inserting ( [13.14| ) into (|13.13|) and taking into account only the first term in (|13.16|) 
we find the Fermi energy of a gas of fermions in curved space-time 



2/3 



'3nv^Y 



1 1 R 
24 



13.17) 



or 



£FiR) 



1 1 R 
24 



(13.18) 



where 



^(0) 



Gir n 



2 

2^\ 3 



\2m 



(13.19) 



is the Fermi energy at the Euchdean space. 

Taking into account the second term in ( |13.16D we get the expression for chemical 
potential 



TT 



/i(T, R) = EpiR) {!- — 



2 / T 



12 \ep{R) 



+ 



1 R 



TV 



T 



24 m2 12 Vel?) 
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or 



KT,R)=^^^'\T) + ^^/P+ ... (13.20) 
24 m"^ 



which describes the exphcit dependence of the chemical potential of the fermionic gas 
on temperature T and the curvature R of space-time. 

Conclusion 

The thermodynamical potentials of quantum bose and fermi gases (as local ther- 
modynamical objects in a curved space-time) were rewritten in terms of the finite tem- 
perature Green's functions of bosons and fermions by means of the local momentum 
space formalism. The phenomenon of Bose condensation was studied and the criti- 
cal temperature of condensation as a functional of curvature was found ( |12.25| ). The 



non-thermal character of Bose and Fermi distribution functions (|12.13|) and ( |13.9| ) was 



shown. The dependence of Fermi energy (|13.18|) and the chemical potential ( |13.2CI| ) of 



a fermi gas on the curvature of space-time was computed for low temperatures. The 
dependence of chemical potential of bose gas on the curvature of space-time was an- 
alyzed. It was found that the temperature is a local thermodynamical characteristic 
of thermal systems in external gravitational fields. 
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PART II 



INTERACTING FIELDS 
AT FINITE TEMPERATURE 

Introduction 

In Part I of this work we studied ideal quantum systems (systems in which parti- 
cles don't interact with each other, but only with external fields) at definite tempera- 
tures. We described statistical properties and thermodynamical behavior of quantum 
ensembles of these particles in an arbitrary curved space-time. However can develop 
a formalism of finite temperature field theory for applications to the systems of inter- 
acting particles (bosons or fermions) and study properties of such systems in external 
gravitational fields. 

Part II is devoted to studying thermal interacting quantum systems in gravita- 
tional fields. In this part we study renormalizability of the finite temperature self- 
interacting scalar \(p^ model, and consider the phenomenon of non-equality between 
inertial and gravitational masses of the boson in perturbative regime at finite tem- 
perature. 
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In Part II the following topics are developed: In chapter XIV A</?^ model at two- 
loop perturbative regime is considered. The concepts of renormalizations and all 
necessary counterterms are described. In chapter XV a complete renormalization 
procedure for finite temperature model X(p^ in two loop approximation of the per- 
turbative scheme is developed. The Green's function of a boson in a heat bath is 
computed. 

In chapter XVI the finite temperature Hamiltonian of a boson is constructed. The 
phenomenon of non-equality between inertial and gravitational masses of a boson in 
non-relativistic aproximation at high temperature in the heat bath is described. 
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Chapter 14 



TWO-LOOP 



RENORMALIZATIONS IN 



A</>* MODEL 

We will start with the problem of renormalization procedure for the self-interacting 
scalar model in two-loop approximation of perturbation theory. In our calculations 
we will use the method of counterterms. 

Let the Lagrangian of the self-interacting Ay?^ model be 

L = l{dv?-lmW-^v' (14.1) 



We may assume that (|14.1| ) is 



L = Lo + Lr (14.2) 
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where the Lagrangian 



Lo = \ - \mW (14.3) 



describes propagation of free particles, and the Lagrangian 



1 - ,2,„2 ,A 



Lj = --6mV - -f^' (14.4) 



describes interaction. 

Let the coefficient Sm"^ be the difference of the form 



c- 2 2 2 

dm = — 



and constant Ab is expressed as 

AB=/i'-"(Afi + 5A) (14.5) 

We determine and tjir as bare and renormahzable boson masses and A^ and \r 
as bare and renormahzable constants of interaction. 

It is easy to see from (|14.3 ) and ( 14.4 ) that the free propagator of a scalar field is 



G = i/ {p^ — mjj 
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and its vertex is 



4—n 




Feynman rules for this model have a standard form ||ltzykson fc Zuber 1980| . 



The Feynman diagrams of the counterterms may be found from the Lagrangian 
of interaction (|14.4|) . 

The two point counterterm diagram is 



and the vertex counterterm is 

Taking into account these counterterms one can construct contributions of the 
order Ab to Feynman propagator G 




a) b) c) 
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Fig. II-l One loop and counterterm contributions to the self energy of the boson. 
Propagator Fig. la) has a standard form. 

The self-energy diagram of the first order to for Fig. 1-b) may be constructed 
from the above mentioned graphs in the form 



Fig.lb = 




1 (-zAR/i^-") 

2 (27r)" 



k'^ — m\ 



1 {Xrh 



4-n\ 



2 (27r)" J ^ '"k^-ml 



(14.6) 



The factor (1/2) is the symmetry factor and is a parameter with dimension of mass. 
This parameter is used to absorb the dimension of the coupling constant. 



Using the expression [[Ryder 1985|| for the ra-dimensional integral 



(Tk 



(/c^ + 2kq — m^)° 



— 1) Z7r2 



,r(«-f) 1 



r (a) (g2 + ^2^ 



;i4.7) 



we find 



Fig.lb 
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2 (27rr(^2)i-t' 2 J 



\ — — 1 



r 1 

2 167r2 Utt/x^ / V 2 



n 



From the expression for gamma function 



r 1 



n 



2J n-4 



+ 7-1 



14.9) 



finally get the divergent part of the self-energy diagram 



Fig. la 




iXr IV? 

IGvr^ 72 — 4 



+ -^i? X finite term + O iXji 



;i4.io) 



The renormalization of the first order of \r may be done with the equation 



: — - + {-i5n?) - ^ 



IGvr^ n — A 



+ 



;i4.11) 



Let us express term 5im? in the form of series ||Collins 1974 



oo oo U \1 
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m 



2 I buXn bi2X'^ b22X% 



R 



+ 



+ 



(n-4) (n-4) (n - 4) 



+ .. 



(14.12) 



where the coefficients bi,j are the numbers. 



In order for vertex T^'^^ in one loop approximation to be finite 



irP) ^p^-rnl+(sm' + ^ fir^Ue (14.13) 

\ 16tt^ n — Aj 



assume 



5m 



2 Aijml 1 



167r2 n-4 



(14.14) 



then 



167r2 



;i4.15) 



This result completes the one loop calculations. 
Now we will find the vertex corrections for F*^^). 
There are four graphs of the order \\: 




+ 




+ 




a) 



b) 



c) 



+ 




d) 



Fig. II-2 Feynman diagrams contributing to the vertex correction F^^^. 



174 



The vertex counterterm Fig. 2a) was introduced before. 

The loop contribution (Fig. 2b)) of the order A|. to the vertex function is 




1 {-iXrH 



4— n\2 



2 (27r)2« J [k'^-m?) {p + kf - m 



2 ,,8-n 



2 (27r)2" 



(/c^ — m^) (1 — x) + + k)^ — m 



X 



2 ,,8-n 



2 (27r) 



2n 



[A;^ + 2pkx — (m^ — p^x)]" 



(14.16) 



Taking into account ( |14.7D we get 




2(27r)2" V 2 



— p'^x{l — x)] 2 



1 

IGtt^ n — 4 



+ A|j X finite term 



(14.17) 



Here we used the representation for F (2 — ?t,/2) in the form 
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f n\ 2 



(14.18) 



Two other loops have the same divergent contributions, therefore F*^^-* vertex structure 
can be expressed in the form 



,2 ,,4-n 



p(4) -x ,A-n 3 A^/i „-x\,.4-n 

J- — — ^A/j/i 



167r2 (n - 4) 



— i5X^ 



r(4) = -i\ -4-" 



3 A^/i 



2 ,,4-n 



167r2 (n - 4) 



+ 5A 



From the equation ( p.4.19| ) we find, that 



14.20) 



Putting 



AB = /i"-"(A^ + 5A) 



4— n 



00 00 \J 



(14.21) 



we find coefficient 012: 



ai2 



167r2 



14.22) 
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Therefore 



Xb = H 



4—n 



(14.23) 



Further we will consider two-loop contributions to two-point vertex F*^^^ 




+ 




+ 




+ 




a) 



b) 



c) 



d) 



Fig. II-3 Counterterms and the loop contributions of the order X% to the self 
energy. 

To find counterterms to the order 0{X^), we will make loop calculation of all these 
contributions. 

Diagram 3(a) gives 



Fig.3a — 




1 (-zA y -") 

2 (27r) 



1 Xnijf' " / 2\ 
272^ ) 



^^^r (2 - 1) 

(m2)'-5 
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(14.24) 



Inserting ( |14.18| ) into ( |14.24| ) we get 



Fig.3a 




+ 



(167r2)2 (n - 4)2 2(167r2)2 (n - 4) 



finite terms (14.25) 



Diagram 3(b) gives 



Fig.3b 




1 {—i6\fi 



4-n\ 



2 (27r)" J k^-m^ 



(Tk- 



l^A/i^-" . . Hxr(i-t 



2 (27r)' 



(m2)^~^ 



(14.26) 



Inserting (|]A9|) into ( pA2^ ) we find 



Fig.Sb 
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+T^7#^Z^(7-1) (14.27) 



(167r2)2 (n - 4)2 2(167r2)2 (n - 4) 

Diagram 3(c) gives: 




Fig.Sc — 



j dl'k j (Tl- 



4 (27r)''^ i J {k'' - rn?f {P - rn?) 
I A|/i^-2" /• 1 f 1 



4 (27r 



4 (27r)'" J ^ '"{k^-m^fJ ^ - m?) 



(/. HxT (2-1)1 _._„,r(i-f 



\ 1 '''t I O * » I '''' / 1 IT' 



4(167r2)2 I m2 / V 2/ V 2 



tXl ml iXl ml 2] + finite terms (14.28) 



(167r2)2 - Ay {32ttY {n - 4) 
Diagram 3(d) is 



Fig.Sd = 
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6 (27r)^" J ^ '"J ^ \k^-m^){P-m^)[{p + k + l)-m^] 



(Tl 



6 (27r)^" J J - w?) {P - w?) [{p + k + l)- w?] 

The function r(3 — n) may be written as 



r (3 - n) = + 7-1 (14.30) 

n - 4 



As the result, the expression for 3(d) will be 



Fig.Sd 




,2 ^2 / ^2 \ 



(167r2)' V^TT/^V (n-4)^ 



' f^_!^ + (^_l)^A (14.31) 



(167r2)' (n - 4) [12 2 
The complete two loop calculations give us the following expression for F^^^ 



i 167r2(n-4) 
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or 



(167r2)2 (n - 4)2 ^ 2(167:^ {n - 4) ' 



+ 



(167r2)2 (n - 4)2 2(167r2)2 (n - 4) 



7 



(167r2)2 _ 4)2 (327r2)2 (n - 4) 



iXWR ( m-R' y-' 1 
(167r2)2 1^ 47r/i2 j (n - 4)^ 



.A| 1 /p2 



(167r2)' (n - 4) [12 2 



(14.32) 



•■n(2) 2 2 , c 2 , ^R'ITIr 

iT^ ' — p — m + dm + 



167r2(n - 4) 



2 A2jm2^ ^ 1 



X2 
■^R 



(167r2)2 (n - 4)2 (167r2)2 (n - 4) 



2 2 

12 2 



(14.33) 



To make r(2) finite in the second order of A^, we will put 



(n-Ay (^-4) (n-Ay 



b22 + 0{X%)} (14.34) 



Then, combining terms in the proper way, we get 



m 



R 



(n-4) 



^11 



167r2 



+ 
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1 + 



1 



\2 



12(167r2)2 (n-4) 



1 + 



1 



X2 



12(167r2)2 (n-4) 



-1 



X n + 



X2 



(n - 4)2 



2(167r 



2^2 



>'12 



+ 



(n 



4) 



(167r 



2^2 



-'22 



(14.35) 



It follows from ( |14.35| ) that coefficient bu is exactly ( |14.15| ). 

Two other coefficients are found from the suggestion that T^"^^ is analytic at n = 4. 
In the result we get that 



-"22 



2(167r2)2 



(14.36) 



and bi2 is the solution of the equation 



\2 

■^R 



(n - 4)2 



Jl2 



12(167r 



2^2 



(14.37) 



It gives 



J12 



12(167r2)2 



(14.38) 



After these calculations we will have 



1 + 



X2 
■^R 



12(167r2)2 (n-4) 



2 ™2 



X n — m 



R 



(14.39) 



or 
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(14.40) 



where (wave function) renormalization constant Z is 



Z = 1 + 



1 



\2 



12(167r2)2 (n-4) 



(14.41) 



From the calculations of this section we found that the bare mass rriB and coupling 
constant for the second order of perturbation theory are 



(n-4) 



1 5 A/ 



167r2 12(167r 



2^2 



+ 



2A| 



(167r2)2 (n-4)2 +^ 



14.42) 



and 



\b = /i"-" A^ - 



\2 



;i67r2) (n-4)""^V^^ 



14.43) 



The expressions ( [14.42| ) and ( 14.43 ) connect non-renormalizable and renormalizable 
parameters of the model, and the model is renormalized in two loop approximation 



of the perturbative regime. 
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Chapter 15 

GREEN'S FUNCTION OF 
BOSON 

IN FINITE TEMPERATURE REGIME 

The aim of this chapter is to construct renormahzable Green's function for a boson 
in a heat bath with a definite temperature. For this purpose we will use the real time 
representation for the finite temperature propagator that will let us obtain necessary 
results in a natural and elegant way. 

We will repeat calculations for the contributions in self energy of the boson in one 
and two loop approximations based on the scheme, developed in the previous chapter. 

In contrast to chapter XIV we will consider that all internal lines of Feynman 
graphs are the finite temperature propagators of the form 
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D{k) = Do{k) + Dp{k) 



i 27r(5(P - m2) 



^2 _ ^2 g, 



/3|fcO| _ I 



(15.1) 



where /3 ^ = T is the temperature. 

Finite temperature calculations don't change Feynmann rules for loop calculations 



[ [Morley 19781 , || Bernard 1974 



Fig lb) gives 



Fig.lb 




(T^O) 



1 {-iXRfi' 



4—n\ 



2 {2ny 



(TkD^k) 



<rkDp{k) 



or in a more compact form 



Fig.lb 




(r = o) 



IHAh/^ 
2 (27r)" 



(TkDpik). 



(15.2) 



The counterterm Fig. 3a) is 
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Fig.Sa = 



or 




Fig.Sa — 

(r = o) 



(27r) 

The contribution of the counterterm Fig. 3b) is 



>^RSm'J -^^Dp{k)Do{k). (15.3) 




(T^O) 



1 (-i5A/x^-" 



2 (27r) 



1 {—iSXn 

2 (2^ 
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and 



FigM 




(r = o) 



Two loop contribution Fig. 3c) may be written in the form 



Fig.Sc — 




(T^O) 



1 {-iXriX 



A-n\2 



4 (27r) 



2n 



<rk / (TqiDoik) + D0{k)f{D,{q) + Dp{q)) 



1 {-zXrIi 



4-n\2 



4 (27r) 



2n 



I d'^k I dy [Dl{k)Do{q) + Dl{k)Df,{q) + 



+2Do{k)Di3{k)Do{q) + 2Do{k)Df,{k)Df,{q) + Dl{k)Do{q) + Dl{k)Df,{q) 



Then 
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Fig.Sc — 

(r = o) 



2 ^4^ r ,,4-n 



Finally the contribution Fig. 3d) will be 




Fig.3d — 

(T^O) 

1 (-iAH/i4-")2 



6 (27r) 



2n 



X 



y rqD^{q)Do{q - p - k) + J (rqDp{q)Do{k - p - q) 
I (rqDp{q)Do{q + p + k) + j (rqDp{q)Do{q - k + p) 



^ J ' " 
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+ I d''qDo{q)Do{q-p-k) + J (FqD^{q)D^{q+p + k) 



or 



FigM = 




+ 



(T = 0) 



2 (27r) 



2n 



j (rkDp{k)x 



X 



(rqDp{q)Do{q + p + k) + j rqDo{q)Do{k+p + q) 



So, we get 



FigM = 




(r = o) 



2 J (27r)4' 



(27r)^ 



,4— n 



^ / tI^^/^^^) ^ ^ / d''qDo(q)Do(k + q + p) 



(27r)'^ 



Now we can find counterterms. 
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Assuming the sum of the finite temperature contributions of the equations ( |15.3| ) 
and ( |15.5| ) zero, 



\2 1,4— n „ I „ 



2 (27r 



(27r) 



we find the expression for 5m? in the form 



5m 



2 



2 (27r) 



(15.7) 



The divergent part of this counterterm will be 



167r2 (n - 4) ' 



(15.8) 



The following counterterm 5\ may be found by the summation of the finite temper- 
ature contributions of the equations (0), ( |15.5| ) and ( |15.6| ) 



-^Mk) L5\ + ^ 
(27r)4 ^\ 2 



n— 4 



(27r)" 



+A 



n-4 



(2vr) 



rqDo{q)Do{k + q + p) 



For zero external momentum p the divergent part of the 5X will be determined by the 
divergent part of the integral 



5X 



n — 4 



(15.9) 
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and the divergent contribution will be 



5A.,„ = -. (15.10) 



The temperature counterterms ( |15.8[ ) and ( p.5.10| ) have the same structure as the 



counterterms which annihilate the divergent parts of zero temperature loop contribu- 
tions. 

It is easy to see that at T = the sum of the loop contribution and the counterterm 
gives zero: 




4-n^ 



2(27r) 



(TkD^ik) - W = 0. 



It leads us to the equation (|15.7|) . 

Loop contributions and counterterm at T = in r*^"^) (Fig.II-2) gives the equation 



2 ,,4-n 



2 (27r; 



— / d''kDo{k)Do{k + p)- iSX = 0, 



which coincides with ( [15.9| ). 

From the mentioned above analysis we can conclude that the following loops have 



the same divergent structure: 
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T = 



T ^ 



and 





+ 




+ 




T ^ 



The Green's function D' (p) of the boson is the object which takes into account 
virtual processes of creation and anihilation of the additional particles when this 
boson moves through the vacuum. 

The graph composing D' (p) may be divided into two distinct and unique classes 
of proper and improper graphs (Fig.II-^ 



^The proper graphs cannot be divided into two disjoint parts by the removal of a single line, 



whereas the improper ones can be disjoint Bjorken fc Drell 1965 
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d'{p) = (>-o + e— ^ ^-o + e— ^ ^ o q + ■ ■ ■ 

Dip) D{p) Dip) Dip) Dip)Dip) Dip) 

Fig.II-4 Green's function D' (p) of the boson as sum of proper self-energy insertions. 
In accordance with Fig.II-4 the Green's function -D'(p) is obtained by the summing 
of the series: 



D' ip) ^ Dip) + Dip) (^]Dip) 



+ Dip) (^]Dip)Dip) (^]Dip) + 



^^\ + iY.ip)Dip) p^-w?-T.ip) ^^^'"^^^ 



In this equation E(p) is the sum of all two point improper graphs. All divergent 
improper graphs and their counterterm graphs (Fig. 11-1,11-3) may be divided into 
two parts of the first and the second orders with respect to \r. We will define them 
as El and E2 self-energy graphs. 

One can write finite contributions in S2 in the form 
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,,4-n 



(27r)^ 



J d^kDo{k)Do{k + p)\ (15.12) 

J p=0 



Let us introduce functions 



and 



Then we get 



II ^.m^) (15.13) 



= / ^^/3(g)^o(g+p). (15.14) 



E2 = A|F^/o(0) + (15.15) 
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where 



Gp = {Fp,{Ip + Io))f^nite (15.16) 



The first equation in ( p.5.16| ) is a scalar product of the form 



(27r)4 ' J (27r) 



and the following one is 



{F„ I,)s^nue = I ^,D,{k) I / ^^D,{q)D,{k + q)\ (15.18) 



(27r)4 ' \i (27r 



finite 



The contribution of the first order in Si may be found from ( |15.2|) . 
This contribution is 



(2vr)^ 

The temperature contribution in the boson's mass0 will have the following form 



^For finite temperature quantum electrodynamics the non-equality of fermionic masses Smp/m ~ 
a{T/m)'^ is described by a similar equation |Peressutti & Skagerstam 1982| 
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= ml + XnFp + \lFpIp{Q) + -fXlGp (15.20) 
As the result the Green's function will be 

Thus we have computed the finite temperature Green's function of a boson in two-loop 
approximation in the form of Feynman propagator with finite temperature dependent 
mass parameter ( |15.2(j| ). 
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Chapter 16 

THERMAL PROPERTIES OF 
BOSON 



In chapter XV we showed that the model is renormahzable in each order of the 
perturbative regime, and found the finite temperature propagator of a boson in a heat 
bath in two loop approximation. We also got the expression for the finite temperature 
mass of a boson. These results may help us to get an effective Hamiltonian of the 
particle and to study its finite temperature behavior in gravitational fields. 



16.1 Effective Hamiltonian of the boson 
in non-relativistic approximation 

After renormalization the pole of boson propagator ( |15.21| ) may be written as 
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E 



2 



(16.1) 



We can rewrite the equation (|16.1|) in non-relativistic approximation in the following 
form 



E = mR[l + XRf{(3mR) + o{Xl)Y 



1 + 



ml{l + XRf{Pmn) + o{Xl)} 



rriR + ^XnTUnfi/Smn) + £^(^1 + ^A^/(/3m^)) ' + o(A^) (16.2) 



Here the function f{y) is connected with the function Fp[y) ( [L5.13| ) in the following 



way 



I r (Pk 1 



2 J (27r)3 £ (e^^ - 1) 



The function /(y) has an asymptotic form (for (y <^ 1)) ||Dolan fc Jackiw 1974|| : 



/(l/) = 777-^ I dx- ^ 



2 J """" e^'s^ - 1 
1 



1 1 



2%2 Svr?/ 



(16.3) 



which is very useful for the analysis of the high temperature behavior of the model. 
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16.2 Inertial and gravitational masses of a boson 

For our following calculations we will consider that the quantum system interacts 
with the gravitational field which is described by the metric 



(16.4) 



where $ is a gravitational potential. 

The second term in ( 16. 4| ) describes a small correction to the Minkowski metric 
which is connected with the presence of the gravitational field 
Landau & Lifshitz 11 , IIMisner et al. 19731 . 



In order to write the Hamiltonian of the boson in the presence of gravitational field 
one can consider that temperature T changes according to Tolmen's law flblmen 1934 | 



T 



To 



1 + $' 



(16.5) 



where Tq is the temperature with $ = 0. 

The finite temperature Hamiltonian with precision to the first leading term of the 



series (|16.3|) will be 



T2 



48 (1 + ^ymn^ 



niR + 



X 



R 



48 m/j(l + $) 



7^2 

° +mR<t> + ... (16.6) 



The last term of the equation ( |16.6|) describes energy of boson's interaction with 
gravitational field. 
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Let us rewrite the Hamiltonian ( p. 6 .61 ) as 



The acceleration of the boson in a gravitational field may be found from the quantum 
mechanical relation 



mH,f\] = -[l-^^JVi, (16.8) 



\rT, 



and the mass ratio will be 



!!k = 1 _ ^li 

rrii 12 rn? 



l--^^ (16.9) 

R 



SO the inertial and gravitational masses of the boson in the heat bath are seen to be 
unequal. 

One can estimate the value of Smg/rrii for some gravitational source. Let the 
source of gravitational field be the Sun ( 1.989 x 10^^ kg) then the relation (|16.9|) 
for the combined boson (Cooper pair with mass rrib = IMev) in the heat bath with 
temperature 300i^' gives the following corrections for non-equality between masses 



Srug Xr T, 



rrii 12 



A X 10"^^ (16.10) 



or 
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10-21 < ^ < 10-17 (16.11) 
rrii 



for the range of the couphng constant 10^^ < A < 10^^. 

From the analysis we made in this chapter one may conclude that thermal inter- 
action of the bosons in a gravitational field causes non-equality between inertial and 
gravitational masses. Non-thermal systems do not demonstrate such properties. 

The calculations for non-equality between inertial and gravitational mass of elec- 
tron were made for thermal quantum electrodynamics by Donoghue 
[ ponoghue et al. 1984|] . His result for massive fermions has the same functional struc- 



ture as the equation (|16.9|) 
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PART III 



NON-LINEAR MODELS IN 
TOPOLOGY NON-TRIVIAL 
SPACE- TIME 

Introduction 

The usual method of generating spontaneous symmetry breaking in quantum field 
theory is to introduce multiplets of vector or scalar fields which develop a nonva- 
nishing vacuum expectation values ||Goldstone 196 1||. This mechanism is not nec- 



essary. The general features of spontaneous symmetry breaking are independent 
of whether the Goldstone or Higgs particle is associated with an elementary field 
or with a composite field. However it is possible to develop a dynamical theory 
of elementary particles in which the origin of the spontaneous symmetry break- 
ing is dynamical ||Nambu fc Jona-Lasinio 196 If . The finite mass appears in analogy 



with the phenomenon of superconductivity. This part III deals with the research of 
Bender et al.l977 l, [ Tamvakis fc Guralnik |, [ Kawati fc Miyata 1981 in the field of 



dynamical symmetry breaking and phase transitions. Models with dynamical mass 
generation are interesting for the following reasons: 
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First: Initially massless fermi fields expose 75 invariance, which is violated by 
the dynamical mass generation. Thus the initial symmetry will be broken, and this 
broken symmetry may appear in experimental measurements [[laylor 1976 



Second: Dynamical mass of particles in such models is a function of a constant 
of interaction of primary fields. Perhaps this will indicate how to solve the puzzle of 
the origin of particle mass. 

Third: Higgs mechanism introduces into the theory additional parameters which 
are connected with Higgs- Goldstone fields. These are additional difficulties of the 



model. IIKaku 1988 



Distinctive and important characteristics of the models with dynamic symmetry 
violation are the introduction of bound states of particles, and calculations in non- 
perturbative regime. Estimations of the condensate may be obtained from the sum 
rule lll^'loratos et al. 198411, or from Monte-Carlo method IIBinder 19791, or for the 



semiclassical instanton solutions of the Yang-Mills equations [|Shyriak 1983 



To develop further the theory one should take into account also the topology of 
space-time ||Toms 1980|| , ||lsham 1978|| , fFord 1980[] . One should consider 3-D and 
4-D spinor models with non-linearity of the type {ipipy. The non-linear Gross- 
Neveu model in 3-D space-time in non-trivial topology is studied in chapter XVII. 
3-D Heisenberg-Ivanenko model in non-Euclidean space-time and 4-D Heisenberg- 
Ivanenko model in Riemann space-time are studied in chapter XVIII. The influence 
of topology and curvature on the generation of dynamical mass of fermions and the 
effects of violation of symmetry in these models are also studied in chapter XVIII. 
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Chapter 17 



NON-PERTURBATIVE 



EFFECTS 



IN GROSS-NEVEU MODEL 

Quantum Field Theory may be essentially symplified in the limit of high interanal 
symmetries such as 0(N), SU(N) and so on. In some appropriate cases field models 
can be solved strictly in the limit of large fiavor numbers N [[Gross fc Neveu 1974 1. 



Asymptotically free Gross-Neveu model without dimensional parameters in the La- 
grangian is one of such models. This model is renormalizable in 3-D dimensions. 
The solution of this model shows that the phenomenon of dimensional transmutation 
Coleman fc Weinberg 197^ has place and there is a gap in mass spectrum of the 



model. Studying this model one will find the connection between topological char- 
acteristics of space time and the behavior of the solution of the mass gap equation 
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and also the behavior of dinamical mass of the model. N-flavor Gross-Neveu model 
is described by the Lagrangian 



9 



;i7.i) 



This Lagrangian is invariant under the discrete transformations: 



(17.2) 



The generating functional of the model 



Z[ri,f]] = I DipDipexp i ( (Tx {iip dip + {1/2) g'^{ipipf + f]ip + ipr]) 



may be rewritten in the form: 



Z[r], 77, cr] = j DipDipDa 



X exp 



i / d''x Uipdip - g{^p'(p)o- - (l/2)cr^ + r]ilj + ipr] 



;i7.3) 



Here we introduced the new field a and used a useful relation 



J Da exp —i{l/2) {a,a) — i (^g'iljilj,aj oc exp {i/2)g'^ (^ip^ 



(17.4) 
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Since 



(17.5) 



we must put m = for symmetry of the model. 

Therefore a new Lagrangian may be written in the form 



1 

2? 



17.6) 



and the symmetry of (|17.6|) is 



ip 75?/;, i) -75^^, (J ^ -a 



(17.7) 



The generating functional of the model after integration over the matter fields will be 



Z[0] = J Dip Dip Da exp i J L„ {x)dx 



I 



Da-Det{id — a) exp[-— ^(cr, a)] 
^9 



:i7. 



Then effective action is written as 



r[a] = —i In Det{id — a) -{a, a 

'^9 



;i7.9) 
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and the effective potential is 



(21)" ^ ' 29! 

The conditions for the energy to be minimal are 



/ |ct=(Tc n / |o-=o-( 



>0 (17.11) 



From these conditions one can get the gap equation for definition of cxc (o"c defines 
the minimum of the effective potential) in the form 



I ^ r (Pk I , , 

- = Trl / 17.12 



where constant A = g'^N . 

17.1 Trivial case. Euclidean space time. 

Let us find the solution of the gap equation (|17.12D. 



Ultraviolet cut-off of the integral gives: 

cPk 1 



/. 



A {2iif + a2 



l/27r / j;^-— = l/47r In— 17.13 
Jo k^ + cr^ 
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Then, after regularization of ( |17.13| ) we get the equation 



1 1 

Let the subtraction point be /i = a, then the renormahzed couphng constant may be 
written as 



1 1 

In^ (17.15) 



To ehminate the parameter n one can use the methods of the renormahzation group. 
The /^-function in one loop approximation is 



/?(A^(^)) = -i(A^(/i)f (17.16) 

TT 



then Gell-Mann Low equation 

dX^ 1 dS, 



(17.17) 



with initial condition A (^o) = determines the behavior of the coupling constant 
with respect to scaling of the momentum: 



where t = ln(^/,^o) 
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The dynamical mass can be found from ( |17.11| ) in the form 



acitriv. 



/iexp( 



TT 



Ail 



TT 



yuexp I — T— I = const. 

AQy 



(17.19) 



17.2 Non-trivial topology of space time. 



In this case the constant of interaction can be written as 



of. 



(17.20) 



where /(o"c, V) is some function which depends on topological parameter L. 
The /3 function is 



/?(A^(/i,L)) = -i(A^(/i,L))2 

TT 



(17.21) 



and the solution of Gell-Mann-Low equation with (^qj -^) = Ao(-f') is expressed by 
the equation: 



\r it,L) 



1 + (Ao(L)/7r)t 



(17.22) 



The dynamical mass is governed by the equation 
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or 



(Te = a.itrtv) exp [-1^^^] (17.24) 



One can see that the dynamical mass depends on function /(o"c, L). 
The exphcit expression of the function f{ac,L) is 



oo ^ 1 

= ±^ / (exp ^x^ + {La,f - l) (17.25) 
vr J ^ /3;2 _|_ (^Lu^Y \ J 







for topologies of the cylinder (+) and the Mobius strip (— ). Therefore non- Euclidean 
structure of space time leads to redefinition of the gap equation (|17.12|) for the Gross- 



Neveu model. That gives us the possibility to define the dependence of the fermionic 
mass on topology of the space time. Dynamical violation of the 75 symmetry oc- 
curs when a is not equal zero. The method developed above is very useful when the 
number of "flavors" of the fundamental fields is big (A^ — > 00). This method does 
not work for a small "flavor" number N . In Quantum Field Theory there is another 
method, based on an analogy with superconductivity. This is a Mean Field Method 
[ [Kawati fc Miyata 1981| |, fZinn-Justin 1989| , which works very well for any number of 



"flavors" of the particles. The idea of the method is based on effective potential^ cal- 



^EfFective potential is the generating functional for (IPI) Green's functions | |Jackiw 1974 , 
Rivers 199C | 
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culations, that allows us to take into consideration the effects of topology |[lbms 1980 
and curvature for self- interacting and gauge models [ Ishikawa 1983|| . 



In this work the Mean Field Method is used in dynamical modeling of the behavior 
of elementary particles and is based on the idea that the masses of compound particles 
(e.g. nucleons) are generated by the self-interaction of some fundamental fermion 
fields through the same mechanism as superconductivity. Here the combined particles 
are treated as the quasi-particles excitations. The Mean Field Method also leads to 
the mass gap equation, and the solution gives the dynamical mass of the particle. In 
the following section we will treat the problem of non-Euclidean space-time structure 
in models with a dynamical mass. 
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Chapter 18 

(iptp)'^ NON-LINEAR SPINOR 
MODELS 



18.1 Dynamical mass and 

symmetry breaking 

In the construction of the unified models of the elementary particles one can 
admit the possibility that the mass of combined particles appears as the result of 
self-interaction of certain fundamental fields, for instance, quarks and leptons from 
preons, or compound fermions in technicolor models || Rizzo 1983|, P?erazawa 198CI|| , 



[ [Ito fc Yasue 1984| , || Far hi fc Susskind 1980|| . Following this idea we can get the gap 



equation. Its solution can predict the dynamical mass of the compound particles. As 
in the previous case of Gross-Neveu model, we will study here the effects of non-trivial 
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topology and, also, geometry of background space-time. 

In this section we will consider the phenomenon of the dynamical generation of 
mass of particles in the application to the models in non-trivial space-time. 

Let us consider Heisenberg-Ivanenko [[Ivanenko 1959|] non-linear spinor model with 



Lagrangian 



2 

L = ^idijj + ^ Ui/jf (18.1) 



where is a massless parameter and parameter /Xq has dimension which is connected 
with the dimension of space time. The symmetry of the Lagrangian ( |18.1| ) is 



V-^TsV", ^^-hn (18.2) 



We can rewrite this Lagrangian in the new form 



= iJidi; - goa [iP^) - (18.3) 



The symmetry of the last one is ( |17.7| ). 
The equation of motion for the a field is 



a = U^) (18.4) 



thus we may assume that the field cr is a collective field. 
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Let us consider that a is a — > cr + a, where 



is the background field and a is quantum fluctuations of the cr-field. 
Then 

2 2 

= ^Pitd - g,a)^ - go (^a^) - - - /ioVa (18.5) 



The last term of ( |18.5| ) may be eliminated because of the redefinition of the sources 



of the quantum field a, and the Feynman graphs will be 



ik - goao (V /^o) 
a 

b) c) 
Fig. III-l Feynman graphs including a collective field. 

Graph la) describes the interaction of fermi field with collective field, 16) is the 
propagator of fermi field, and Ic) is the propagator of collective field. 

In tree approximation with respect to cr-field we can write an effective action: 



^e/zM = (-^2) \riDet{ik - goo) - 




(18.6) 
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The effective potential of tliis model will be 



Veff = ^(Tri) / T^ln {k^ - (goar) + f (18.7) 



(2 

Minimum Veff gives the gap equation 



o"c = {go/^iofaj {gocTc) (18.8) 



or, in another form, 



m = sXofnl (m) (18.9) 

where m is the dynamical mass of fermionic field m = goac, s is the dimension of 7- 
matrices and Aq = ((^q/z^o)^ 

Now one can solve the equation ( p.8.9 ) for different space time topologies. 



18.2 Model with topologies 

Ri X Ri X Si and Ri x Mobius strip 

The solution of this gap equation is connected with the calculation of the function 
I(m) of the equation (|18.9|) . 



Let us consider two types of topologies: 
1) Ri X Ri X Si with iplx, y, 0) = iplx, y, L) 

and 
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2) Ri X Mohius strip with ip{x,y, 0) = 
We can find for 3-D space time tliat 



-ij{x,y,L) (Fig.III- 2) 



I{m) = A - F(±){L,m 



18.10) 



wfiere 



F(±)(L,m) = — ln(l;e-^'") 



18.11) 



witfi (+) for Ri X Ri X Si topology and (-) for Ri x Mohius strip topology. 
The gap equation will be 



m = mXoA (l ^— ln(l 7 e"^"*)) 

V TiLA J 



;i8.12) 



The analysis of this expression can be made by the theory of bifurcations 
[Hale fc Kocak 1991|| . For this purpose let us write the gap equation ( |18.12| ) in the 
form 



m = f(±){m, A) 



18.13) 



where the functions f(±){rn,X) for topologies (+)and (— ) are 



/(±)(m, A) = mA ( 1 - ^ ln(l ; e"^™) 



18.14) 
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The equation ( |18.13| ) has stable m=0 solutions, if /m(0, A) < 1. 



For /m(0, A) > 1 the equation ( |18.13| ) has no stable trivial solutions. 
One can see that there is a stable m=0 solution for (-) topology with the critical 
parameter 



n2 Ao n2 
Lc = , / . - — r 18.15 

TT AnA - 1 TtA ^ ^ 



for AqA = a > 1 

The gap equation for topology (+) has no stable trivial solutions. 

The dynamical mass in this case is a smooth function with respect to parameter 

L: 



m(+) = In - /(A) exp (^-^) ) (18.16) 



for AqA = a < 1 

The solution for topology Ri x Mohius strip is 



m(„) = -yln('exp('^ln2) - 1) (18.17) 



The restoration of the symmetry takes place when L = L^. In this case the condensate 
function of the bound state equals zero 



= ^ < >= (18.18) 
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Now we can consider the models with the more comphcated space-time structures of 
Klein bottle and Torus topologies. 
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18.3 Torus topology Mg) = i?i x x 5i 

and topology M3) = Ri x Klein bottle 

We can introduce topologies -Ri x 5*1 x 5*1 and -Ri x Mohius strip as an identification 
of space points for wave functions Fig. III-3: 

^{x, 0, 0) = L) and ij{x, 0, 0) = -ipix, L, L) (18.19) 

For simplicity consider that these topologies have only one parameter L = L' . The 
gap equations for the topologies will be 



m = m 



\s\l + v^m/(A2L)(7r)-3/2[^ -j=K_,i,{2Lmn) 

I n=l V 



1/4 

+ K_i/2{2Lm^/n^ + P) (n^ + ]\ (18-20) 

n,fc=l 



and 



1 



m = mXs \ 1 + Jm/{K^L){TT)-^'^[Y, ^=K^i/2{2Lmn 



n=l V 2n 



1 //I 

+ (-l)'i^-i/2(2LmVn2 + P) (n^ + ^ ] ^ (18.21) 

n,A;=l 



The solutions of (|18.20|) and ( |18.21| ) give dynamical masses for these topologies 
Kulikov fc Pronin 1989| 1. 



The results of this paragraph show that the application of the Mean Field Method 
to non-linear models gives non-renormalized solutions, though we can obtain some 
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information about the influence of topology on dynamical mass behavior. The evap- 
oration of condensate and restoration of chiral symmentry proceed in different ways 
and are ruled by topology. There are topologies in which these phenomena do not 
take place. We believe that these results are important in the bag models, because 
the energy of bag is dependent on non-perturbative effects and boundary conditions 
IGuidry 1991| . 

18.4 Non-linear spinor (^^)^ model in 

Riemann space-time at finite temperature. 

In this section we will treat the problem of dynamical mass generation of the non- 
linear spinor model in 4-D Riemann space-time at finite temperature. We will get 
the finite temperature effective potential and find out information about the influence 
of curvature of the background gravitational field and temperature on the value of 
dynamical fermionic mass. 

Let the total Lagrangian of the model be 

L = Lg -\- Ljn 

where the gravitational Lagrangian i^ 

"'^For clear understanding of the problem, we study the fermionic system in a weak gravitational 
field. 



(18.22) 
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and the Lagrangian of matter field Lm is 



gi n2 



= ^tr^^^ + ^ UV- J (18.24) 



is a covariant derivative. 



The first term of ( |18.24 ) describes kinetics of the fermi field and its interaction 



with the gravitational field, and the second one describes the interaction of the fields. 
The effective potential of the model is written from the action 



r,ff[a] = -'-\nDet[V' + \R-{gocrn - ^{a,a) (18.25) 



After making ultraviolet regularization and renormalizations of the model we get 



Ltot = Lg + Lm = (Lg + L„(oo)) + LmiP) 



Uf^ + LmiP) (18.26) 



The renormalized gravitational constant Gr and the A/j-term are: 



1 r 1 r 1 .o 1 



87rGp 87rGn lovr^ 16n^ 
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where A is a cut-off parameter, and m is a dynamical mass. 

In the calculations ( p.8.26| ) and (|18.27|) there was used the cut-off method of reg 



ularization of divergent integrals Pogoliubov fc Shirkov 1980 



The effective potential can be written from (|18.25|) in the form 



where 



and 



'18.2^ 



2m 



2 CO 



E 



/ q In 1 + 



PJ (27r) 



;i8.29) 



(18.30) 



The coefficients aj{R) are 



ao = 1, ao = —R, 



;i8.31) 
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The solution of the gap equation 



^K//M|.=.. = (18.32) 



gives, in high temperature approximation ( PCXII.33D , the expression for dynamical 



mass without redefinition of the coupling constant and temperature 



m\R, T) = {goa^Y = A/ X + b ■ + C ■ R + ... (18.33) 



where the constants 



A = 327rV3.84, B = A/24:, C = 0.7/12, A = (go/f^o) 



are numerical positive coefficients. As we can see from the equation ( |18.33|) , the 



effective dynamical mass of fernion is a positive function for any temperature and 
curvature. 
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PART IV 



TOPOLOGICALLY MASSIVE 
GAUGE THEORIES 

Introduction 

3-D non-Abelian gauge models describing systems of point particles carrying non- 
Abelian charge have been under investigation for over two decades. These models (so 
called topologically massive models) have a number of interesting features: 

1) For vector fields, these models possess single, parity-violating, massive, spin 1 
excitations, in contrast to single, massless, spin exitations in the Maxwell theory, 
and to a pair of spin 1 degrees of freedom in gauge non-invariant models with a mass 
Pinegar 1982 | 

2) For second-rank tensor fields, describing gravity, the topological model leads to 
a single, parity- violating, spin 2 particle, whereas a conventional (gauge non-invariant) 
mass term gives rise to a spin 2 doublet. Furthermore, the topological term is of 
third- derivative order, yet the single propagating mode is governed by the Klein- 
Gordon equation. Einstein gravity, which is trivial and without propagation in three 
dimensions, becomes a dynamical theory with propagating particles. 
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3) Particles interacting via the Abelian Cherm-Simons term (CS-term) acquire 
anomalous spin and fractional statistics. They are called anyons [Wilczek 1977 1. 



Anyons play a role in the fractional quantum Hall effect [ |Arovas et al. 1985[| , [Lee & 



Fisher 1989] and perhaps also in high temperature superconductivity [Lauglin 1988]. 

This is not a complete list of interesting properties of gauge theories in an odd 
number of dimensions. An important part of these models is the Chern-Simons ac- 
tion (CS-term). Some interesting aspects of quantum field models arising from the 
topology of odd- dimensional manifolds are discussed in chapter XIX. The origin of 
the CS-term of vector type that is induced by gauge interaction of 3-D fermions is 
studied in chapter XX. In this chapter the influence of topology and temperature 
effects are also considered. The origin of the induced gravitational CS-term at finite 
temperature is considered in chapter XXI. 
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Chapter 19 



INTRODUCTION 

TO TOPOLOGICAL FIELD MODELS 

As an introduction to odd dimensional topological field models let us consider their 
origin and topological significance. 

Topological aspect of the model 

Prom gauge-invariant fields in even dimensions we may construct gauge invariant 
Pontryagin densities: 

P(2) = -{l/2n)€'"'*F'"' 
P(4) = -(l/167r^)tr*F'^^F^, (19.1) 
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whose integrals over the even dimensional space are invariants that measure the topo- 
logical content of the model. 

These gauge invariant objects can also be written as total derivatives of gauge 
invariant quantities 

Pn = d^Xi: (19.2) 
The two and four-dimensional expressions are 

= (l/27r)e'^^A, 

X^ = {l/2n)e''''''Hr{A^Ff,^ - {2/3)A^Af,A^) (19.3) 

The Chern-Simons (CS) secondary characteristic class is obtained by integrating one 
component of X^ over the (n — 1) dimensional space which does not include that 
component. 

Therefore the 3-D action Scs is proportional to 

Scs ^ j dx^dx^dx^Xl (19.4) 

A topological massive term (we will name it CS term) can be added to the fundamental 
action for a gauge fields, but unlike the ways in which gauge fields are usually given 
a mass, no gauge symmetry is broken by its introduction. 
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Quantum aspects of 3-D field theory 

Let us consider non-Abelian quantum model with topological mass term. The 
Lagrangian of this model is 

L = Lo + Lcs + Lgauge (19-5) 

where Lq is the usual action for non-Abelian gauge field 

Lo = -(l/2)tr(F^,F'^'^) (19.6) 

with 

F^,, = d^A, - d,A^ + g[A^, A,] (19.7) 

Lcs is CS term 

Lcs = -ime^''nT{A^d,A, - {2/3)A^A,A,) (19.8) 
and Lgauge includes the gauge-fixing term 

Lgauge = {& ^v'') {d^ V,'') + Q f abc{d ^.T,'') T] (19.9) 

We introduce SU{N) gauge group here with matrix notation: = ■^It"'^ where r" 
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are anti-Hermitian matrices in the fundamental representation: 



[t\ = P^'^t", tr(TV) = -(1/2)5' 



ab 



(19.10) 



and /"^"^ are the structure constants of SU{N). 

The theory is defined in three space-time dimensions with Euchdean signature 
(+ + +). The couphng of the CS term is imaginary in Euchdean space-time and 
real in Minkowski space-time. 

Properties of the model: 

For an odd number of dimensions, the operation of parity, P, can be defined as 
refiection in all axes: 



p 



—X' 



(19.11) 



The usual gauge Lagrangian is even under parity. 




(19.12) 



but the CS term is odd 




(19.13) 
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Under gauge transformation 



A^^^-'{{l/g)d, + A^}^l (19.14) 



The Lagrangian Lq is invariant but Lcs is not: 



jd^xLcs^ Jd^xLcs 
+ iim/g) j d'^xe^'''Pd^ii[{d^n)Vt-^Af] + Sti"^ {m / g'^)iuj (19.15) 



where: 



(1/24) j d^xe^"'Pii\VL-\d^^)VL-\d^n)VL-\dpVL)] (19.16) 



The set of gauge transformations is divided into global gauge rotations, d^Vt = 0, and 
all others, for which we assume that Q{x) 1 as ^ oo. Integrating over global 
gauge rotations requires the system to have a total color charge equal to zero. In 
this case, v4^(x) falls off faster than l/\x\ as x^ —>■ oo, and the second term on the 



right-hand side of (|19.15|) , which is a surface integral, vanishes. 



The last term in ( |19.15|) does not vanish in general. The u of (|19.16| ) is a winding 



number, which labels the homotopy class of f2(x) [[Nakahara 1992 | 

Changing variables A A'^ , where A'-' is a gauge transformation of A, implies 
that the vacuum average of the value Q is 
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< Q >= exp \i87i'^{m/g'^)u{U)\ < Q > 



(19.17) 



This invariance gives us a quantization condition for the dimensionless ratio 



Therefore, for the theory to be invariant under certain large gauge transformations 
(for a non-Abehan gauge group), which are not continiously deformable to the iden- 
tity, the ratio of the CS mass m and the gauge couphng g"^ must be quantized 
I Deser et al. 1982 . 

Further we can study the problem of massive excitations. 

For spinor electrodynamics in three dimensions we have 





(19.18) 



L — Lg + Lf + Lint 



(19.19) 



where 




(19.20) 




'int 



(19.22) 



231 



The coupling constant e has dimension (mass) 
The equations of motion will be 

(id + eA — m)ip = 
One can introduce the dual field strength tensor in 3-D space-time 

The Bianchi identity follows from ( |19.23| ): 



di,*F^' = 



and the equation ( |19.23| ) may be written in a dual form 



da*Fp — dp* Fa — IJ-FafS — —eap^J^ 



or 



This equation demonstrates that the gauge excitations are massive. 
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3-D gravity. Connection with topology 

The results of the topological part of the introduction gives us a hint how to 
construct the topological term for three dimensional gravity from a four- dimensional 
*RR Pontryagin density. 

*RR = {ll2y-^^R^,p,R^p'"^ = d^X^ (19.29) 



Let us find from (|19.29|) . To do this we will rewrite *RR in the following way 



*RR = {l/2)e^''^''R^,,^R^f,P^ = (l/2)e^'^"^/?^,,,i?, 



ab 



= e^'^^fR^^at {9,^/^ + iv^^'^up,'} (19.30) 
In our calculations we used the following expression for the Ricci connection 

Rfiupa = dai^^ab — d^^aab + '^aa'^lic}, — ^Pa^^acb (19.31) 

where u^j^ab is 3-D spin connection. 

Then the expression for Pontryagin density will be 

*RR = e^^"^i?^,,, a^cu/" + e'^'^^^i?^,,, uja^'^uj^c' (19.32) 
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Let us find these two contributions separately. 
The first contribution to ( p.9.32| ) gives 

- 6^""''a;/3'^'9„ {d^co^ab - d.co^ab + u;''^,u;,cbtot,cu^,b} (19-33) 
= e'^^^^a^li?^,,,^/^) - 2e'^^"^a;/''9„ {co'^^co.cb} (19.34) 
The second contribution to ( |19.32| ) will be 

f R^uab 



2e^'"'''d^ij,abOJo.'''ojp,'' + 2e^'"'''ujl^u,,,uj^'^'up,' (19.35) 



On the other hand the second one to (|19.35D is zero, then we have 



"RR — e^^'^^ da{Rfj,vab^l3°'^) 
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- 2e'^'^«^a;/^9„ {^^.cb} + 2e>'^''Pd^uj,abU^^'^ujp,'> (19.36) 
One can find that 

= 2e>^^^Pd,ul^uJ,''up,^ (19.37) 

and 

2e'^^"^a^u;^„u;«6'=u;^e" = (2/3)e'^^"^a^ (19.38) 
In the result we will have 

*RR = a^e''^"^ {uJ^^^'Ra^ab + (2/3)cu^„cu«,^cu^e"} (19.39) 

Therefore 

^ [uj,-^R^^^, + (2/3)u;,aVb'a;/3c"} (19.40) 

Let parameter ji be equal to zero, then determining e^''"'^ — e^"'^ with u,a,P — 1, 2, 3 
we get the CS action in the form 

Scs ~ / d^xX'' - J d^xe"'''' {uj^'^'Ra/iab + {2/3)u;Juj^b'ooi3c''} (19.41) 
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As we can see from ( |19.41[ ) the CS term is of the third derivative order in contrast to 



the first one as in the vector case ( |19.3| ). 
Nonlinear theory of gravity 

We can construct total gravitational action in the form 

Stot = (l/k^) J d'x^R + (19.42) 

The sign of the Einstein part is opposite to the conventional one in four dimensions. 
The Einstein part of the action has coefficient with the dimension of mass, while 
the topological part has a dimensionless coefficient (yU has a dimension of mass). 

Now we can find some interesting properties from this action. 

Varying (|19.42|) with respect to the metric, we get field equations 



Q^"" = C"" + (1 /fi)C'"' = (19.43) 



where the second rank Weyl tensor C'^'^ is 



c^" = {i/VaV'-^D^K (19.44) 



and 



RaP = Rap-{ll'^)9apR, R=K (19.45) 
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The components of the Einstein tensor G""^ are 



G"^ = i?"^ - (l/2)c/"^i?, (19.46) 



and the components of the Riemann tensor R^p-^s are 



R^p^s — d&^°'i3-i — d.-^^'^iis + ^""n-y^^ps — ns^^p-y (19.47) 



From ( |19.43| ) we get first order form for field equations 



i^^/"(/i)i?A. = (19.48) 



where K^^^" is operator of the form 



K,.^'' = {6i6t - {l/2)g,,9'n + ^6/^(5^5,^ - img'^g^^) (19.49) 



Operator K^J'''{^) may be multiphed by K^J'''{—^) to yield a second-order equation 
for Ricci tensor. 



From ( |19.48|) we find 



ir„/'^(-/i)ir^/'^(^)i?A. = (19.50) 



that is 
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(19.51) 



This exhibits the massive character of the excitations. 

After this short review of topological field models we will consider 3-D fermionic 
models interacting with vector and with tensor fields. These models may give effective 
induced topological action of the CS type. 
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Chapter 20 



INDUCED CHERN-SIMONS 



MASS TERM 



IN TOPOLOGICALLY NON-TRIVIAL SPACE-TIME 



In gauge theories with fermions, the a topological mass term is induced by fermionic 



interactions [| Redlich 1984|] . If the topology of space-time is not trivial, topological 
parameters will appear in the effective action and in the CS Lagrangian. In this 
chapter we show that the CS term depends on topological parameters in such a way 
that the topological gauge invariance of the total action can be maintained. 

Let us consider the case of a massive fermion field interacting with external gauge 
field in 3-D space-time. 
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20.1 Euclidean space-time. Trivial topology 

The Euclidean action of this quantum system is given by 



S= d^xip{x){idx + m + gA{x))ip{x) 



(20.1) 



The corresponding generating functional is 



Z^l D^{x)D^^[x)exv [- / Sxi.{x){^^. + m + gA[xmx) 



(20.2) 



Here A^(a;) = A^(a;)T", and are generators of the gauge transformations. We will 
choose to work with two-component Dirac spinors. 
Euclidean Dirac matrices have the algebra: 



(20.3) 



This equation is satisfied by the matrices 



(20.4) 



where o", are the Pauli matrices. 



Furthermore, they satisfy 



240 



liup 



pXu) 



The integration in (|20.2|) over the fermionic fields gives 



Z[A^,] = Det{id + m + gA)) = exp{-Seff) 



Therefore we can write the expression for the effective action 



S^jf = - In Det{id + m + gA)) = -Tr \n{id + m + gA) 



Now one can spht ( p0.7| ) in two parts 



5'e// = — Trln(i9x + m + gA{x)) 



Trln^j-Trln 1 + gSfA 



where 



Sf 



1 id + m 



id + m d"^ + m? 



is a fermionic propagator. 
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The first term of the equation ( |20.8| ) is a divergent gauge-independent contribu- 
tion, and the second one may be expanded in a power series [|Babu et al. 198'^: 



S, 



-Tr In 



OO 1 

l + gSfA\=TTY^-[gSfA] 



n=l 



n 



(20.10) 



Let us consider the term which is quadratic in A^{x). The corresponding action is 
given by 



c{2) _ 1 2rp 



id + m -i id + m -> 
-A— -A 



Q2 _|_ ^2 Q2 _j_ ^2 



(20.11) 



The terms with two and four 7 matrices contribute to the wave function renormahza- 
tion of the gauge boson and will not be taken into consideration in further calculations. 
The other two terms are linear with respect to mass m and give 



i9 J 1 J 1 J i9 J 
:A— -A + — -A— -A 



Q2 _|_ ^2 Q2 _|_ ^2 Q2 _|_ ^2 Q2 _j_ ^2 



img 



trz)7M7i^7ptr 



-A 



:An H — ~Aiit::: ~Ar. 



Q2 _|_ ^2 f}2 _|_ j^2 (}2 _|_ ^2 f}2 _j_ ^2 



(20.12) 



where Tr = tr£)tr is the trace of both the 7 matrices and the vector field, and also 
includes the integrations. 

To write tr in the form tr [f{d)g{xy\ one can use the relation 



1 W + J^ + E^4# + ... ,20,3) 



d"^ + m? d"^ + m? (S^ + m^)^ (^^2 _|_ ,^2^3 
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Then 



img e^^p 



(^2 + m2)2 



\x > tr{{d,A^)Ap} (20.14) 



In another form 



S^^^f = -img^ e^,yp / d?x 



d?k 1 
(27r)3 + m2)2 



(20.15) 



The additional contribution to the CS term comes from the terms which are cubic in 



74^'s. We can rewrite them as 



id + m id + m id + m i 
:A— -A-z -A 



Q2 _|_ ^2 _|_ Q2 _|_ j^2 



(20.16) 



Noticing that the terms having even number of 7 matrices do not contribute to the 
CS term, we have: 



S% = -\g^mTr 



:A— -A— -A 



52 -|- 777,2 Q2 _|_ j^2 Q2 _|_ j^2 



^ ^ 1 2 9 

Q2 -1-77^2 Q2 _|_ j^2 ^2 _|_ j^2 
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+ — -A— -A— -A 



-m 



-A 



+ m' 



-A 



Q2 _j_ ^2 



(20.17) 



Using the identity (|2al3|) find 



-.(3) _ 

'e// - o 



3 r 



X tr 



d^dp AAA 9^0\ AAA ^v9\ 

-AyA\A„ + -j7^ — ; ^T^AjyApAf^ + :T7^Ai,Ar,An 



(9^ + m'^y 



(9^ + m^) 



(9^ + m 



1 



^AuAyAr 



Tracing with respect to 7 matrices, find 



(20.18) 



S. 



(3) _ 2mgf^ 



(27r)3 (P + m 



2^2 



tr y rf^xe'^'^M^A^Ap (20.19) 



(20.15 


) with ( 


20.19 


) we get 



term in the form 



= -img'' 



(27r)3 (F + m 



2^2 



tr / (Txe 



dpAyAp — —igApA,yAp 



(20.20) 



or, after integration over the momenta 



qCS _ ig"^ m I ,3 



d/xAiyAp — —igApAjyAp 



(20.21) 



Now we can apply the method we have developed for a gauge model in 3-D non-trivial 



space-time. 
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20.2 Non-trivial topology 

a) Model with space-time topology S = i?^^) x mobius strip. 

To introduce topology S = R^'^^ x with parameter we rewrite the expression 
( |20.15 ) in the following way. Momentum integration in ( P0.15 ) will be integration with 



respect to two dimensions {{t, x) i?^^)),and the third integration will be transformed 
in the sum, because the fermionic propagator is antisymmetric with respect to the 
selected axis {y S^). 

We can consider that k"^ = k'^+ool with Un = (27r/C)(n+l/2) and n = 0, ±1, ±2, .... 
Integrating in 2-D momentum space we find that 



d^k 1 r d^k 1 ^ 1 



^ tanhr-^] (20.22) 



8n\m 



where 



d^k ^ '' 



(2vr)3 Cn^J (2vr)2 
In the process of calculations we used the useful equation 



y 
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Then the expression ( |20.20| ) in this topology will be 



Mcs) . , r d'k 1 ^ 



tr j e^,p \d,ApA, - '^gA.A^A,^ (20.24) 
and the induced CS term at non-trivial 3-D space-time is written as 

Sef%) = --£j^\tanh |M| tr I d^xe,,,{d,A,Ap - hgA^A^A^) (20.25) 
The relation between CS terms is the function of the form 



CS I 

(I I I f ^ r I / I 

(20.26) 



b) Model with space-time topology S = i?*^^) x 5*^. 

For this space-time topology the propagator of the fermionic field will be periodic 
at the boundary points of interval [0,C], that leads to the modification of equation 
( |20.21 ) with respect to new frequencies Un = 21171/ Q with n = 0, ±1, ±2, .... 



Taking into account the summation formula 



E 71^ = ^<^oth{-Ky) (20.27) 



we find, that 
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d^k 1 



d^k 1 g 



(2vr)3 (P + m2)2 y (2vr)2 ( ^f-^ + ^2 + ^2y 



57r m 



(20.28) 



Using this equation and the equation ( p0.19|) , get the induced CS term in the new 
topology 



SeffiO = -^S^o^^ tr / d'xe,,,{d,A,A, - '^zgA.A^A^) (20.29) 



iTT \m\ 



The rotio of the CS terms will be the function of the topological parameter (: 



(20.30) 



These results show that the relations (|20.26|) and (|20.30|) are smooth functions of the 
topological parameter (. 
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Chapter 21 



GRAVITATIONAL 
CHERN-SIMONS 



MASS TERM 
AT FINITE TEMPERATURE 



In this chapter we will consider fermions interacting with an external gravitational 
field at finite temperature. In the previous chapter XX we learned that the interaction 
of fermions with external gauge bosons leads to the induced action of the Chern- 
Simons type. In the same way the interaction of fermions with external gravitational 



fields may leads to the effective gravitational Chern-Simons term || Redlich 1984 



Alvarez-Gaume et al. 1985| . Now we will develop a formalism of calculations of the 



effective gravitational CS action from the action for massive fermions interacting with 
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an external gravitational field. 



21.1 Induced gravitational Chern-Simons 
mass term 

Let us introduce the action for massive fermions connected with an external grav- 
itational field as: 



S = j (fx^i){x){iD + m)i){x) (21.1) 

where D = = 7^(9^ + uOab^j.cr"'^) with 7^(x) = h^{x)'y"-. 

The function uJabfi is the local Lorentz connection and o""^ = (1/8) [7"^, 7^] is the 
commutator of 7 matrices. 

7-matrices with Latin indices are constructed from Pauli matrices 

= ia" ,a''{a = 1,2,3) (21.2) 

and obey the relations: 

{7^ 7'} = -25'^^ 7,7b = -5ab - eabcl" (21.3) 

tr(7'^7'') = -2i(5°^ tr(7'^7%") = 2e'^''" (21.4) 



We can find the effective action by integrating ( |21.1| ) with respect to the fermionic 
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field from tlie equation 



exp[-Seff] = J D'ip{x)ip{x) exp[-S] = Det{-iD - m) (21.5) 



To find the one-loop effective action we will rewrite ( pi .51 ) in the following form 
Ojima 1989 : 



Seff = —In Det ^ 



In Det- 



iD + iuj — m) 
—iD — m){—iD — iuj — m) 



{—iD + iuj — m)(—iD — iuj — m) 
= - In Det{DD - Duo + imuo + m^) + logDet[{D - u)^] + (21.6) 

where we put uj = = ■-y^a°'''ujab^i- 

Now we can use the momentum space formalism we developed in part I of this 
work for calculation of the first determinant of the expression ( ^1.6] ). We can simplify 



the problem by noticing that the second term of ( |21.6| ) does not give contributions 
to the induced CS term. 

For the calculation of the first term let us introduce normal coordinates and write 
the metric at the origin (x) of these coordinates. We will use variable y = x — x for 
the definition of the point of manifold in tangent space. The tangent space Y will be 
"flat" with metric gab = —i]ab = diag{—l, —1, —1). In the origin of these coordinates 
the tetrad functions and the metric will be the series: 
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Kix) = - (i/6)i?^,.i/V + - 



(21.7) 



The covariant derivative in a tangent space is written as 



D = i,{x)D^{x) 



rK{x)-{i/Q)rR'.a.y'y'' + ..- 



'A 



+ cu„b^(x>'^M (21.9) 



The effective action can be represented in the normal coordinates as 



S. 



eff 



J d^x^/gSeff 



X 



hm, / rfV V^5e//(x') = hm / d'x^Seffix, y) (21.10) 



In proper time formalism InDet is 



Seff[oj\ = In Det{DD — Du + imuj + 



hm^ j d^x < x\ J — trexp j-i (j)''D^ ~^ T ~ ^'^ ~^ ^^'^ ~^ "^^) } ^ 



In accordance with ( |21.1CI| ) and ( |21.11|) the density Seff{x) may be written in the 
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momentum-space form 



/ds f d^Iv 
— J j^—j^ exp[ iky ] exp[-i(m^ - A^A^)s] x 



tr exp < —i 



Df'D^ + — - Diu] + {iA^D'' + iD^Af" - iAcu + imuj) 



s\ (21.12) 



where lirn A^(?/) = fc^. 

Notice that if we hke to obtain the contribution of finite terms as (m — > 0), it 
is sufficient to extract the terms which are proportional to with a factor m or 
proportional to with a factor (or A;^). Thus we get the term in Seff[uj], which 
contains the induced CS term. 

The second and the third orders of the expansion ( pi.l2| ) give 



/d^kt f ds 

— ---exp[ iky ] / — exp[-z(m^ - A^,A^)s](zm)x 

^ ' 



X tr 



2! 



{ujDuj + Dujuj) s — 



3!^ ' 



(21.13) 



where A{y) = 7''A^(t/). 

In the limit {y 0) the equation ( pl.l3|) is 



/d^ k f ds 
(2^ J — expH(m2 + A;^)s]x 



252 



X tr 



UJUJUJ s 



— — {kukujuj + kuLukuj + ujkuikui + m ujujui)s 



(21.14) 



Noticing, that 



ii(kujkuLiuj + kujuokuj + uokookuj + m uoujuj) = {k + m )tTUJUJUJ (21.15) 



rewrite ( pi.l4| ) as 



Seff{x) = {im) 



iiuoduj + -trcjcjcu 



Sk 
(27r)3 



rfs(s)exp[-is(m^ + r)] (21.16) 



Applying the equation for the r(z/) function 



Tiu) 



\m 



2u 



(21.17) 



and computing tr of 7 matrices, find 



I m 



64:TT \m\ 



J d^xe'^'^P {d,uj\,uj% + '^uj\,uj\,uj\p^ (21.18) 



This is the expression for gravitational CS term |Pjima 1989|| which we get usin^ 
momentum-space method. 
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21.2 Induced gravitational Chern-Simons 
mass term at finite temperature 

The calculations of finite temperature gravitational action of CS type are based 
on the computation of In Det (|21.5| ). In Part I we developed the formalism of such 



calculations with the help of the momentum space methods. For these calculations we 



will replace the integration procedure over three momenta in the expression ( |21.16| ) 
by the summation over Matsubara frequencies and integration over the two remaining 
momenta. 

Let us introduce temperature in the tangent 3-D space to the curved manifold as 



(part 1,(113)): 



/ (0^^^''^''^^'''^^L£^/ ^2nu^n,klkls) (21.19) 

where = T is temperature, and = (27r//3)(n + 1/2). 

To introduce temperature in the model, rewrite integral over (s) in ( pi.l6| ) in the 
form 

IWfl'' -'^Pl-l*^ + - - / (0^FT^ (21.20) 



Using (|21.19|) and the summation formula ( p0.23|) , find the finite temperature CS 
term 

7 771 771 \ r / 2 \ 

Seff[u^] = J d'x'''"' [9^^'auu^\ + -u;\,u',,u\^j (21.21) 
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Comparing (|21.18|) and (|21.21D we get the following result 



Seff[T = Q] 2 



(21.22) 



The expression (|21.21|) shows that the structure of the induced CS term at finite 
temperature is exactly the same as at zero temperature, and the functional relations 
between the CS terms in tensor (|21.22| ) and vector types of interactions ( |20.26| ) are 



the same. 
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APPENDIX 
SUMS OF BESSEL FUNCTIONS 

I. Type I sums of the modified Bessel functions 

The integral representations for series of modified Bessel functions may be calcu- 
lated from the following summation formula: 



E 



n + l 



-1 



TT 



tanh(7r?/). 



(XXII. 1^ 



In proper time representation the left side of (|XXII.1|) can be written in the form 



E 



2 I n + 1 



°? oo 

/ da exp (^—ay'^^ ^ exp 



a [n + 



(XXII.2) 



Taking into account the equation 



exp |— a (n — 2;)^! = ^ i — j exp ( n^ — 2mzn\, (XXII. 3) 
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we can write ( pCX11.2| ) as 



E 



2 I n + 1 



CO 1 
'TT \ 2 



a 



-n 



a 



- + 2 V (-1)" Ida ( -\ ' exp i-ay^ - -J 



(XXIL4) 



The right side of ( pCXll.l| ) is 



TT vr 2-7? 

-tanh(7r?/) = TT^^TV^- 

y y y [f^v + 1) 



(XXII.5) 



Therefore, we get from ( pCXIL4|) and ( [XXI 1. 5| ) the following useful equation 



z(e^ + 11 



1 oo °? 1 



n=l 







(XXII.6) 



Moreover, we may consider that = Qaix"^) is the function of variable x & E? with 
a parameter a, namely = + a? . 

Integrating with respect to x one can get 



X 



(2vr)^ 



\/ x"^ + a? {exp {\/ x"^ + a^^ + 1^ 



1 CX) „ / 

^ n=l n V 



-72 



(XXII. 7) 



The modified Bessel function may be written as 
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da ■ a'^ ■ exp —7a 



(XXII.8) 



Then ( PCX11.7|) will be 



l/a\^(-ir 



(XXII.9) 



ra=l 



Scaling x and a with a parameter /5 ) = (3{k,m) write ( pCX11.9| ) in the form 



(27r)' £ (exp + 1) (37T^ 



n 



(XXII. 10) 



where e = y + w?. 



Differentiating the equation (PCXII.7|) with respect to parameter a we get 



d^x ( d 



(27r)' V 



\/ x"^ + a? (^exp {\/ + a?^ + 1 



-| 00 

.— E(-ir^o(an), 

n=l 



(XXII. 11) 



or, in new variables, 



f d^k ( d \ 1 ^00 



(XXII.12) 
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Integrating in ( pCXll.7] ) with respect to parameter (a) and using the equation 



a2 + 1 



2 In (l + exp + 



(XXII. 13) 



we find with new variables the following equation 



1 r d^k A (-1 



ln(l + exp(-/3£:)) 



^ ^-J-K2{f3mn). (XXII. 14) 



2(/3vr)2„tl n2 



High temperature asymptotes [jSm <^ 1) of the equations (PCXII.9| ), (pCXII.ll| ) and 
(PCX1I.14| ) are 



2m2 - (-1)- 

2^ —^K2{l3mn) 



7 H ^ + -m^ [In- h 7 (XXII. 15) 



^ E - + (m (^) + 7 - (XXI1.16) 



£lE(-ir^o(/3mn) = iflnf^V7V 



27r 



?i=i 



(XXII.17) 
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II. Type II sums of the modified Bessel functions 



Let us start with the sum: 



E (y 



' n^l = — cothfvry). 

y 



(XXII. 18) 



The proper time representation the left side of ( pCXlI.l^ ) will be 



oo 1 I" °° 

^ (y'^ + n'^y = daexp(^-ay^) ^ exp(-an2). (XXII. 19) 



Taking into account the equation (pCXII.3| ) and putting {z = 0) we find 



exp |— an^l = ^ ( — j exp | n 



1 



(XXII.20) 



Then ( PCXILIQ]) may be written in the form 



^ ' y + n 



n=—oo 



-1 



da { — ] exp I —ay'^ 

aj \ a J 



77 ^-^ f 1 I I 2 ^ 2 

— I" 2 2^ I day — ] exp ( —ay n 



y 



n=l 







a 



a 



(XXII.21) 



The right side of (PCX11.18| ) is 



- coth(7ry) = - + — ^- — - 
y y y (e^'^f - 1) 



(XXII.22) 
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Therefore, we get from ( pCX11.21"D and ( pCX11.22[ ) the following useful equation 



Z TT 



1 \ ^ r /7r\2 
— ; = — 7: > / da — exp I —a — n" 



(XXII.23) 



Let = ga{x'^) be a function of variable x & B? with a parameter a of the form 



After integration with respect to x we get 



d^x 
(2^ 



Vx^ + (^exp (^\/ x'^ + 



a2 - 1 



(XXIL24) 



- ^ / daa ^ exp - 
2 n=i r{ V 



-a—— n 



(XXII.25) 



The modified Bessel function is written as 



/ 



da ■ d"-^ ■ exp f -7a - ^ j = 2 V iT^ ( 2^57 



(XXII.26) 



and ( PCXII.251) will be 



d^ 



X 



^/xP~+~a? (exp (^\/x^~+~a?^ — 1 



1 / a \ ^ 1 



Kiian) 



(XXII.27) 
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Scaling x and (a) with a parameter /3 as (x, a) = f3{k, m) write ( pCX11.27D in the form 



(27r)' e{exp{(3e)-l) (3ti^ 



m °° 1 



(XXII.28) 



where e = v k'^ + m? . 



Differentiating the equation (PCXII.25D with respect to parameter a we find 



\/ x"^ + (^exp (^yj x"^ + a?^ — f^^ 



^" n=l 



(XXII.29) 



or, in new variables, 



d^k ( d 



(27r)' \dm'^) e (exp {(3e) - 1) 47r2 



Y,Ko{pmn). (XXII.30) 



Integrating (|XXII.25|) with respect to parameter (a) and using the equation 



oo 

J do? Vx^ + (^exp (^Vx^ + 



2 -f 



2 In (l - exp + 



(XXII.31^ 



we find with new variables the following equation 



,2 oo 



^E^AV/Smn). (XXII.32) 
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High temperature asymptotes (/5m ^ 1) of the equations ( pCX11.27D , ( |XX11.29D and 
(^Xig) are 



2m^ °° 1 



S 4 

m m 
90p^ 24/32 ^ Unp ^ 64 



167r2 2 ' 



(XXII.33) 



m m 



24/?2 87r/3 327r2 



In— ^-1 + 27 



(XXII.34) 



and 



1 



+ 



1 



167rm/3 327r2 



In 



+ 27 



(XXII.35) 
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GRAPHICS 



gs/2{z,R) 



93/2{l,R) 
93/2{l,R) 




Fig. I-l Graphical expression of the function g3/2{z, R) 
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gs/2{z,R) + X^no 



^3/2(1, i?) + A^rio 

2.612 
^3/2(1, i?) + A^rio 




R<0 
R = 
R>0 



Z\ Zq Z2 



Fig. 1-2 Graphical solution for bosons. 
Solution of the equation (|12.21|) for different curvatures and fixed temperature and 
density 2:1 for _R < 0, zq for _R = and Z2 for i? > 
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T' Tr T'l 



T 1 r 



i?<0 i? = Oi?>0 



Fig.I- 3 Chemical potential /leffiR) as a functional of a curvature of space-time. 
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f3/2{z,R) 




^1 ^0 ^2 



Fig. 1-4 Graphical solution for fermions. 
Solution of the equation ( |13.13|) for different curvatures: 2;i for i? > 0, zq for i? = 

and Z2 for i? < 
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Ri X Ri X Si 



Ri X Mobius strip 




Fig. III-l Topologies of cylinder and Mobius strip in {y, z) spaces. 
Identification of the points y, 0) = y, L) and y, 0) = y, L) 
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Fig. III-2 Topologies of torus and Klein bottle in {y, z) spaces. 
Identification of the points '4^{x, , 0) = '4>{x, L, L') and -0(3;, 0, 0) = —i^{x, L, V) 
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